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PREFACE 


This  paper  represents  the  culmination  of  one  phase  of  the  research 
done  by  Anne  Kanle  during  the  twelve  years  she  was  on  the  staff  of  The 
Rand  Corporation.  During  a part  of  this  time  she  was  enrolled  in  the 
doctoral  program  at  UCLA  and  this  paper  is  her  doctoral  thesis.  During 
the  period  when  this  work  was  done  there  was  a close  collaboration  be- 
tween Rand  scientists  and  the  staff  at  UCLA. 

Although  this  paper  is  primarily  concerned  with  the  solution  of  a 
difficult  physical  problem,  many  of  the  results  were  used  for  the  solu- 
tion of  operational  problems  presented  to  Rand.  This  paper  is  being 
issued  in  order  to  make  these  results  more  widely  available. 
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SUMMARY 


The  problem  of  radiative  transfer  in  a plane-parallel,  perfectly 
scattering  atmosphere  is  described.  Chandrasekhar's  solution  applicable 
to  atmospheres  of  small  and  moderate  optical  thickness  is  outlined.  His 
solution  reduces  the  problem  to  that  of  determining  the  X-,  Y-,  K- , and 
L-functions,  the  scattering  functions.  Mullikin  has  extended  this  method 
of  solution  to  atmospheres  of  large  optical  thickness. 

Sekera  and  Kahle  have  used  Mullikin's  method  of  solution  for  calcu- 
lating the  emergent  radiation  from  plane-parallel  Rayleigh-scattering 
atmospheres  of  large  optical  thickness.  Their  numerical  results  are 
reproduced  here  in  the  Appendix,  as  tables  of  scattering  functions.  The 
numerical  method  for  determining  the  intensity  and  polarization  of  the 
radiation  emerging  from  the  top  and  bottom  of  atmospheres  is  given,  and 
suggestions  for  additional  uses  of  the  tables  are  made.  Finally,  a few 
examples  of  representative  calculations  are  presented. 
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INTRODUCTION 

One  of  the  basic  problems  of  radiative  transfer  in  a planetary 
atmosphere  is  that  of  determining  the  intensity  and  polarization  of  the 
radiation  emerging  from  the  top  and  bottom  of  the  atmosphere.  The  logi- 
cal first  step  in  this  problem  is  the  study  of  scattering  in  a pure 
molecular,  or  Rayleigh  scattering  atmosphere.  One  of  the  most  success- 
ful and  widely  used  methods  is  that  of  Chandrasekhar.  His  solution  for 
plane-paral lei  atmospheres  involves  solutions  of  nonlinear  integral 
equations  for  his  X-  and  Y-functions  by  successive  iterations,  start- 
ing with  the  solutions  for  small  optical  thickness.  This  method  was 
applied  by  several  workers.  However,  when  they  attempted  to  extend  the 
solutions  to  atmospheres  of  large  optical  thickness  (x  > 1 ) the  solu- 
tions failed,  due  to  computational  instabilities. 

Recognizing  that  in  dealing  with  planetary  exploration  we  might 
have  to  deal  with  atmospheres  of  large  optical  thickness,  an  attempt 
was  made  at  the  Rand  Corporation  to  extend  these  calculations  to  large 
optical  thickness  by  using  the  solution  for  x -*-  <»  as  a starting  value. 

However,  the  calculations  failed  to  converge,  and  it  became  evident  that 
the  solutions  were  oscillatory.  Mullikin  made  an  extensive  mathemati- 
cal study  of  the  nonuniqueness  problem  and,  by  going  back  to  the  origi- 
nal radiative  transfer  equations,  was  able  to  select  the  correct  solu- 
tion [Mullikin,  1962a, b,c,d;  1964a, b].  He  devised  a method  of  solution 
which  was  both  unique  and  computationally  desirable  as  its  rate  of  con- 
vergence increased  as  one  looked  at  larger  optical  thickness.  Carlstedt 
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and  Mullikin  (1966)  computed  the  X-  and  Y-functions  for  large  optical 
thickness  for  isotropic  scattering. 


By  appropriate  modification  of  their  method  Sekera  and  Kahle 
(1966)  were  able  to  use  this  same  method  for  Rayleigh-scattering  at- 
mospheres of  large  optical  thickness.  The  tables  of  X-  and  Y-functions 
and  the  related  K-  and  L-functions  they  derived  are  given  here  in  the 
Appendix.  These  constitute  the  complete  solution  to  the  problem  of 
radiative  transfer  in  a plane-parallel  Rayleigh-scattering  atmosphere 
of  any  optical  thickness,  to  as  great  an  accuracy  as  desired.  Kahle 
then  examined  some  aspects  of  the  intensity  of  radiation  emerging  from 
Rayleigh  scattering  atmospheres  of  large  optical  thickness  [Kahle, 
1968a],  and  also  the  global  radiation  fields  for  the  same  problem 
[Kahle,  1968b]. 

In  Chapter  I basic  definitions  are  given,  and  the  problem  of 
radiative  transfer  in  a homogeneous  plane-parallel  Rayleigh-scattering 
atmosphere  is  posed.  The  description  of  the  problem  follows  that  of 
Chandrasekhar  (1950)  and  is  presented  here  as  a review  of  the  concepts 
and  terms  necessary  for  an  understanding  of  the  subsequent  chapters. 

The  solution  of  this  problem  by  Chandrasekhar's  method  employ- 
ing X-  and  Y-functions  is  discussed  in  Chapter  II.  In  Section  A the 
solution  as  developed  by  Chandrasekhar  (1950)  is  examined.  Section  B 
describes  Mull  ikin’ s extension  of  this  solution  to  atmospheres  of 
large  optical  thickness.  In  Section  C the  method  used  by  Sekera  and 
Kahle  (1966)  to  calculate  the  solution  for  a Rayleigh-scattering  at- 
mosphere is  described. 
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In  Chapter  III,  the  tables  of  X-  and  Y-  and  K-  and  L-functions 
given  in  the  Appendix  and  their  use  are  described.  Section  A describes 
the  values  available  in  the  tables.  The  equations  required  to  compute 
the  Stokes  parameters  from  these  tables  are  given  in  Section  B.  Other 
uses  for  the  tables  and  the  accuracy  of  the  tables  are  discussed  in  the 
subsequent  sections. 

In  the  final  chapter  two  applications  of  the  tables  as  solu- 
tions to  the  radiative  transfer  problem  are  briefly  examined  as  an 
illustration  of  their  use. 

In  retrospect,  one  might  well  ask  what  use  can  be  made  of  these 
calculations  now.  It  was  recognized  at  the  time  the  computations  were 
made  that  the  atmospheres  of  Venus  and  Jupiter  were  cloudy.  However, 
it  was  believed  at  that  time  that  there  well  might  be  extensive  niol ocu- 
lar atmospheres  above  the  cloud  tops,  where  these  computations  would  be 
applicable.  Care  would  have  to  be  exercised  in  this  application  since 
the  reflection  from  thick  clouds  would  probably  not  be  perfectly 
Lambertian.  Although  we  now  know  that  sucli  a molecular  atinosphere  does 
not  exist  above  the  Venusian  cloud  tops,  these  calculations  may  still 
be  applicable  to  Jupiter  where  there  appears  to  be  approximately  2 bars 
of  H2  atmosphere  above  the  cloud  tops. 

Another  reason,  at  the  time,  for  making  these  calculations,  was 
the  concept  that  non-molecular  effects  could  be  considered  as  a per- 
turbation on  molecular  scattering.  This  argument  lias  its  limitations, 
however,  because  of  possible  non-linear  interactions  between  the 
molecular  components  of  scattering.  Subsequently  several  authors  have 
attempted  to  solve  the  problem  for  cloudy  atmospheres,  for  both 
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homogeneous  and  non-homogeneous  cases.  In  all  such  studies,  numerical 
methods  of  solution  have  been  required.  There  is  always  a need  for 
testing  such  models  against  a standard  calculation,  the  results  of 
which  can  be  justified  on  analytic  grounds.  It  is  hoped  that  the 
present  tables  will  provide  such  a standard,  since  the  results  of 
these  more  complicated  models  should  reduce  to  the  results  of  the 
Rayleigh-scattering  model  under  the  appropriate  conditions.  This 
should  provide  a necessary  but  by  no  means  sufficient  condition  that 
such  models  should  satisfy. 
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I. 


THE  PROBLEM  OF  RADIATIVE  TRANSFER 


In  this  chapter  the  problem  of  radiative  transfer  in  a plane- 
parallel  Rayleigh-scattering  atmosphere  is  outlined.  Section  A con- 
tains definitions  of  some  of  the  principal  quantities  used  in  radiation 
studies.  In  Section  B the  Stokes  parameters  for  polarized  light  are 
defined.  In  Section  C the  radiative  transfer  equation,  defining  the 
change  in  the  Stokes  parameters  on  interaction  of  the  radiation  field 
with  a Rayleigh-scattering  atmosphere,  is  presented.  Finally,  in  Sec- 
tion D,  the  scattering  and  transmission  matrices  are  defined,  the  quan- 
tities whose  determination  will  constitute  the  solution  of  the  problem. 

A.  BASIC  DEFINITIONS 

The  basic  problem  of  radiative  transfer  is  to  determine  how  a 
radiation  field  is  altered  on  passing  through  and  interacting  with  a 
medium.  If  we  consider  a pencil  of  radiation  of  intensity  I^  (in  a 
frequency  interval  v to  v + dv  ) propagating  through  a medium  in  a 
specified  direction,  then  we  can  consider  the  various  ways  this  inten- 
sity will  be  altered  in  travelling  a distance  ds. 

The  radiation  can  be  weakened  by  two  processes,  true  absorption 
by  the  medium,  and  scattering  of  the  radiation  into  another  direction. 
In  practice,  both  of  these  processes  are  combined  into  a mass  attenua- 
tion  coefficient,  defined  as  < in 

V 


A 
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where  p is  the  mass  density  of  the  medium.  This  attenuation  coeffi- 
cient can  be  due  to  true  absorption  only,  scattering  only,  or  a com- 
bination of  both.  The  case  of  scattering  only  will  be  considered 
first.  If  dm  is  the  mass  of  the  element  of  the  medium  under  consid- 
eration, then  the  rate  of  scattering  of  energy  out  of  the  element  is 

K I dm  dv  dw  (1-2) 

where  duj  is  the  solid  angle  of  the  incident  pencil  of  radiation.  The 
rate  at  which  energy  is  scattered  into  a given  solid  angle  djj'  can  be 
found  if  we  specify  the  phase  function  for  scattering,  p(cos  ■,-)  where 
9 is  the  angle  between  the  incident  radiation  and  d;.. ' . This  rate 

is 

K^I^P(cos  0)  ^ dm  dv  d.o  (1-3) 

The  total  amount  scattered  will  thus  be 


K I dm  d 

V V V 


p(cos  0) 


4n 


(1-4) 


which  shows,  by  comparison  with  Equation  1-2,  that  we  must  normalize 
the  phase  function  such  that 


p(cos  0) 


(1-5) 
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To  account  for  absorption  we  use  the  same  expression  (1-4), 
only  now  do  not  require  that  the  phase  function  be  normalized.  Instead, 


we  have 


- 


(1-6) 


This  defines  uig,  the  albedo  for  single  scattering.  This  quantity  is 
a measure  of  how  much  radiation  has  been  scattered,  while  (1  - ujg)  is 
the  amount  absorbed.  When  Wg  = 1 we  have  perfect  scattering. 

The  phase  function  p(cos  -J)  depends  upon  the  type  of  scatter- 
ing process  involved.  The  simplest,  p(cos  G)  = constant,  is  for  i so- 
tropic  scattering.  The  phase  function  for  Rayleigh  scattering,  of 
great  interest  in  atmospheric  problems  being  appropriate  for  molecular 
scattering  of  visible  radiation  (or  any  scattering  by  dielectric  par- 
ticles which  are  small  compared  to  tfie  wavelength  of  light  both  outside 
and  inside  the  particle  [Van  de  fiulst,  1957]),  is  given  by 


p(cos  fl)  - (1  + cos^O) 


(1-7) 


Rayleigh  scattering  is  perfect  scattering,  i.e., 


yi)(cos  0)  (1  + cos^fi)  ^ - 


= 1 


(1-8) 
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As  with  the  weakening  of  our  pencil  of  radiation  in 


traversing  an  element  of  mass,  there  are  two  similar  processes  enhancing 
the  radiation:  scattering  of  radiation  into  the  direction  of  the  pencil 

and  emission  by  the  mass  element.  Again,  the  two  are  treated  together, 
this  time  with  a single  emission  coefficient  defined  such  that  the 
element  of  mass  dm  emits  into  the  solid  angle  dc.  an  amount  of  radia- 
tion in  the  frequency  range  v to  v + dv  , 

dm  doj  dv  (1*9) 

in  unit  time. 

The  contribution  to  this  emission  due  to  scattering  from  a di- 
rection (0',ii')  into  the  direction  0,<))  will  be  (cf.  Equation  1-3) 

dm  dv  dcj  p(e,c^:0’<})')  1^(0' ,!?>■)  -4-/ ^ 0-10) 

Integrating  over  all  incoming  angles,  the  emission  coefficient  due  to 
scattering  only  is 

j (0 j4>)  = < T” 

•'v'  V 4tt 

(Ml) 


2tt 


I f 


p(0,4>:0.4> 


)I  (0' 

' V 


,4>‘ ) sin  0'  d9'  d^' 
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In  general,  of  course,  there  will  be  true  emission  also,  but  we  will  not 
concern  ourselves  with  that  in  this  paper.  The  source  function  is  de- 
fined as 


J E ^ 

V < 

V 


(1-12) 


For  our  case  of  scattering  only 


TT  2tT 


h f f .<}>')sin  0'  de'  d4)' 


0 0 


(1-13) 


The  basic  equation  of  radiative  transfer  sums  up  the  con- 
tribution of  the  various  processes  which  we  have  just  discussed 


di"  = ■ ^ V 


(1-14) 


, dl 

_^  = I . j 

K^p  ds  V V 


(1-15) 
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The  source  function  is  usually  an  integral  function  of  the  inten- 
sity (see  Equation  (1-13))  so  this  equation  is  usually  an  integro- 
differential  equation. 


It  is  the  solution  of  this  equation,  for  the  goometry  of  a homogeneous 
plane-parallel  Rayleigh-scatteri ng  atmosphere  that  is  the  subject  of 
this  study. 

A homogeneous  plane-parallel  atn.osphere  is  defined  to  be  an  at- 
mosphere which  is  stratified  in  plane  parallel  layers  so  that  the  only 
variation  of  atmospheric  properties  is  in  the  vertical  direction.  This 
variation  can  be  incorporated  into  the  normal  optical  thickness,  t , 
generally  measured  inward  from  the  top  of  the  atmosohere  (considered  to 
extend  to  infinity). 


z 


(1-17) 


The  radiative  transfer  equation  then  becomes 
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(1-18) 


dl^ 


where  u = cos  6 

B.  STOKES  PARAMETERS  FOR  POLARIZED  LIGHT 

In  order  to  describe  completely  the  nature  of  the  radiation 
field  in  the  atmosphere  we  must  look  not  only  at  the  intensity  of  the 
light,  but  also  at  the  state  of  polarization.  The  system  chosen  by 
Chandrasekhar  (1950)  and  his  many  followers  to  describe  the  polarization 
in  the  radiative  transfer  problem  is  that  of  using  the  Stokes  param- 
eters. These  were  first  introduced  by  Stokes  in  1852  [Shurcliff,  1962]. 
These  parameters  are  best  suited  to  both  the  theoretical  and  experimen- 
tal description  of  incoherent,  partially  polarized  light  which  is  the 
natural  occurrence  when  sunlight  interacts  with  planetary  atmosplieres 
[Oeirmendji an,  1969].  All  four  parameters  have  the  same  physical  di- 
mensions, the  parameters  for  coincident  streams  are  additive,  and  they 
are  relatively  easily  measured  experimentally.  The  radiative  transfer 
solutions  of  Chandrasekhar  for  the  intensity  of  the  radiation  can,  in 
general,  be  applied  to  the  Stokes  parameters  by  merely  substituting  the 
appropriate  matrices  for  the  phase  functions  and  other  quantities. 

The  polarization  of  electromagnetic  radiation  is  traditionally 
described  in  terms  of  the  electric  vector.  If  a^^  and  a^  arc  the 
scalar  components  of  the  vector  electric  field  in  two  mutually  perpen- 
dicular directions,  perpendicular  to  the  direction  of  propagation,  then 
the  Stokes  parameters  arc  defined  as 
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■ ‘ h*‘r 

5'  1^/-  »/:■=  'r 

U = {2aj^a^  cos  y} 

V = {2a^j^a^  sin  y } 


(1-19) 


where  y is  difference  of  the  phases  . The  brackets  indi- 

cate a time  average  over  aii  appropriate  time  length.  All  these  Stokes 
parameters  have  the  dimension  of  intensity. 

The  electric  vector  traces  an  ellipse,  and  it  can  be  shown  that 
the  principal  axes  of  the  ellipse  are  in  directions  making  angles  x 
and  X + j with  the  direction  1 where 

tan  2x  = ^ (1-20) 

Also,  the  ratio  of  the  major  and  minor  axes  of  the  ellipse  is  tan  S 
with  the  sign  of  6 determining  the  direction  of  rotation  of  the  elec- 
tric vector  and 


sin  2R  = 


(Q^  + 


(1-21) 


Sometimes  and  are  used  in  place  of  the  first  two  Stokes 
parameters,  I and  Q . 
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For  completely  polarized  light,  I = Q + U + V . For  par- 
2 2 2 2 

tially  polarized  light  I > Q + U + V and  the  percent  of  polari- 
zation 


P 


^ 100(Q^ 


(1-22) 


For  natural,  or  unpolarized,  light  Q = U = V = 0 . For  linearly 

polarized  light,  V = 0. 

If  I = I^,  U,  V),  then  a rotation  of  axes  through  an 

— V 

angle  $ will  subject  I to  a linear  transformation. 


2 2 
cos  ({)  sin  $ 1/2  sin  2;  0 


L M 


2 2 
sin  4-  cos  4>  -1/2  sin  2^  0 


-sin  2({)  sin  24>  cos  24> 


(1-23) 


C.  RAYLEIGH  SCATTERING 

3 2 

The  Rayleigh  scattering  phase  function  p = (1  + cos  >j), 

which  was  mentioned  in  Section  A,  is  appropriate  for  natural  light, 
but  when  considering  polarized  light  we  need  a more  complete  descrip- 
tion of  the  scattering  process.  For  polarized  light  characterized  by 
the  vector 


J 


4 


(1-24) 


r 


T = (i^.  V*  u,  V) 

which  is  incident  on  a single  particle,  the  scattered  light  in  the  di- 
rection 0 will  be  given  by 


dui 


-)  R I duo 


(1-25) 


upon  Rayleigh  scattering  [Chandrasekhar,  1950].  In  this  expression  c 
is  the  scattering  coefficient  per  particle  given  by 


0 = 


2 

a 


(1-26) 


for  Rayleigh  scattering  where  a is  the  particle  polarizability  (to  be 
defined  later)  and  X is  the  wavelength  of  the  incident  liglit.  Also, 

R in  Equation  (1-25)  is  the  phase  matrix  for  Rayleigh  scattering, 
defined  as 


R 


3 

2 


0 

1 

0 

0 


0 

0 

cos  0 


0 


cos  0 


(1-27) 


11 

For  natural  incident  light,  with  I = (^  I,  y I,  0,  0),  R 1 

3 2 

reduces  to  4 (1  cos  0)I,  equivalent  to  the  Rayleiglt  pliase  function 
given  carl ier. 
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So  far,  looking  at  a single  particle,  our  coordinate  system  has 
been  determined  by  the  angle  of  incidence  and  the  scattering  angle.  As 
we  go  to  the  more  general  problem  with  many  particles,  we  must  trans- 
form each  scattering  event  to  a chosen  coordinate  system  applicable  to 

the  problem.  This  can  be  done  by  use  of  the  rotation  transformation 
^ -►  ■* — ► 

L ('{))  given  in  Section  B,  Equation  (1-23).  The  phase  matrix  P will 
be  given  by 

V (e,  4);  O'  , 0' ) =T  (tt  - i2)^R  (cos  6)  I (-  i^ ) (j-28) 

whore  i^  and  i^  are  angles  relating  the  incident  and  scattering 
directions  to  the  chosen  coordinate  system.  It  can  be  shown 
[Chandrasekhar,  1950]  that  "P  can  then  be  written  as  the  sum  of  three 

V - >-  / Q ) 

terms,  an  azimuth  independent  term  P ^ ; a term  dependent  upon 

(<))'  * and  a term  dependent  upon  2(^'  - ^),P 

For  incident  natural  light,  rather  than  incident  arbitrarily 
polarized  light,  the  problem  reduces  somewhat.  The  resulting  scattered 
radiation  will  be  partially  plane-polarized.  Further  scattering  will 
change  the  angle  of  the  plane  of  polarization  and  the  degree  of  polar- 
ization, but  it  will  remain  plane-polarized.  Mathematically,  we  can 

■<  - ► 

say  that  the  matrix  P is  reducible  with  respect  to  V.  Therefore, 
for  sunlight  (natural  light)  incident  on  a Rayleigh-scattering  atmos- 
phere we  need  only  consider  the  three-dimensional  Stokes  vector  and 
three-dimensional  Rayleigh  phase  matrix. 
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The  phase  matrix  for  Rayleigh-scattering  in  a plane-parallel 
atmosphere  with  incident  natural  light  is  thus  given  by 


4pij'  cos  (<J)'  - 0 2y  sin  ($'  - ({)) 


- 2u'  sin  ($'  - (j>)  0 cos  (o'  - o) 


(1-32) 


V m'  .4)' ) 


cos  2 (o'  - 4>)  - cos  2 (o'  - O)  sin  2 (O'  - o)^ 


4|  - y'  cos  2 (o'  - 0)  cos  2 (o'  - O)  - p'  sin  2 (o'  - O) 


yy'^  sin  2 (O'  - O)  P sin  2 (O'  - O)  yp'  cos  2 (o'  - O) 


(1-33) 


The  mass  scattering  coefficient  k is  given  by 


K = - N 
P 


(1-34) 


where  o is  the  scattering  coefficient  per  particle  given  at  the  begin- 
ning of  this  section,  p is  the  mass  density,  and  N is  the  number  of 
particles  per  unit  volume.  For  molecular  scattering  the  polari zabi 1 i ty 


n'-  1 


(1-35) 


where  n is  the  refractive  index  of  the  medium  (cf.  Equation  (I-?6)). 
For  Rayleigh  scattering,  the  mass  scattering  coefficient  is  thus 


< 


8tt^  (n^  - 1)^ 

3 


(1-36) 


This  is  the  quantity  that  enters  into  the  calculation  of  the  optical 
thickness  i 

With  the  phase  function  defined  above,  and  r defined  (cf. 
Equation  1-17)  in  terms  of  given  in  Equation  (1-36),  the  equation  of 
transfer  for  a plane-paral 1 el  Ray  lei gh-scattering  atmosphere  is 


" dr 


1 ( ~ » Vi  ! ) 


+ 1 


1_ 

4ti 


-1 


2m 


P(u,  (i;  ii'  , $'  ) I (t  , p'  , i)' ) dp'  d;  ' 


(1-37) 


THE  SCATTERING  AfJD  TRAllSMISSIQll  MATRICES 

When  considering  the  radiative  transfer  through  a plane-parallel 
atmosphere,  one  of  the  principal  goals  is  to  determine  the  diffuse  ra- 
diation emerging  from  the  top  and  bottom  of  the  atmosphere.  It  is  con- 
venient to  explicitly  separate  the  diffuse  radiation  from  the  atten- 
uated incident  solar  radiation.  If  we  assume  a parallel  beam  of 
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sunlight  of  net  flux  -nF  = F^,  F^)  incident  on  a plane  parallel 
atmosphere  in  the  direction  (-  Uq,  4>q),  then  the  radiative  transfer 
equation  can  be  written 


(x,p,$)  = I (i 


+ 1 2tt 
-1  0 


u',  '+>')  I (t,  p',  $')  dy’  d$' 


1 ”^^^*0-1  >■  -y 

e P(p,  - Pq,  $q)f 


(1-38) 


With  the  direct  solar  radiation  thus  separated  from  the  diffuse 
reflected  and  transmitted  radiation  we  can  define  a scattering  matrix 
S and  a transFiiission  matrix  T as  follows.  Let  the  optical  thickness 
be  measured  from  the  top  of  the  atmosphere;  ■ = 0 at  the  top  and 
T - T.|  at  the  bottom  of  the  atmosphere.  Then  the  reflected  intensity 
is  given  by 


1 ^ • 
^ (0>  4*)  " 'yT!'  S (^1*  4**  4^n)^ 


and  the  transmitted  intensity  by 


I(ti;  - Pt  <J))  = T P>  ‘t'i  Mq* 


0’  ^0' 


(1-40) 


The  solution  of  our  radiative  transfer  problem  is  then  to  find 
+->■  -<— *■ 

S and  T for  a Rayleigh-scattering  atmosphere. 
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II.  SOLUTION  OF  THE  RADIATIVE  TRANSFER  PROBLEM 
BY  USE  OF  X-  AND  Y-  FUNCTIONS 


In  this  chapter  the  method  of  solution  of  the  radiative  trans- 
fer problem  for  a plane-parallel  Rayleigh-scattering  atmosphere  by  use 
of  X-  and  Y-  functions  is  presented.  Section  A describes 
Chandrasekhar's  method  as  he  developed  it.  When  various  attempts  at 
computation  of  the  X-  and  Y-  functions  for  atmospheres  of  large  optical 
depth  failed,  due  to  numerical  instability,  Mullikin  investigated  the 
problem.  He  found  that  the  usual  reduction  of  the  original  radiative 
transfer  problem  to  the  solution  of  equations  for  X-  and  Y-  functions 
resulted  in  the  loss  of  i nfon:iation , causing  nonunique  solutions.  By 
adding  constraints  derived  from  the  original  radiative  transfer  prob- 
lem, he  was  able  to  develop  a method  to  select  the  correct  solution 
from  the  family  of  solutions,  and  thus  extend  the  X-  and  Y-  function 
method  of  solution  to  all  optical  thicknesses.  This  is  outlined  in 
Section  B.  In  the  final  part  of  the  chapter.  Section  C,  the  numerical 
method  of  calculation  of  these  X-  and  Y-  functions,  as  developed  by 
Carlstadt  and  Mullikin  (1966)  and  used  by  Sekera  and  Kahle  (1966),  for 
the  Rayleigh-scattering  problem,  is  outlined. 

A.  CHANDRASEKHAR'S  SOLUTION 

Chandrasekhar  (1950)  has  shov/n  that  the  scattering  and  trans- 
mission  matrices,  S and  T can  be  written  in  terms  of  pairs  of  X- 
and  Y-  functions.  This  separates  the  variables  in  the  problem:  where 
S and  T are  functions  of  (i,  p,  ip,  u',  <(>')  the  X-  and  Y-  functions 
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depend  only  on  (r,  y).  Four  pairs  of  X-  and  Y-  functions  are  re- 
quired for  the  case  of  Rayleigh  scattering.  Isotropic  scattering  re- 
quires only  one  pair.  The  X-  and  Y-  functions  satisfy  a pair  of  simul- 
taneous integral  equations 


1 

Xi(y)  = 1 + y J 

0 


Xi(y)X.(y')  - Y.(y)Y.(y') 
y + y‘ 


dy 


(II-l) 


Y. (y)  = exp  (-  ^)  + y y* 


Yi(u)Xi(y')  - X.(y)Y.(.')  ^ 


M - U 


r ) dy' 


(II-2) 


where,  for  Rayleigh  scattering,  i = 1,  2,  3,  4.  This  notation  differs 
slightly  from  Chandrasekhar , who  used  the  labels  i = 1,  2,  r for 
these  quantities.  The  functions  are  the  characteri Stic 

functions  which  depend  upon  the  type  of  scattering  process.  They  are 
even  polynomials  in  y satisfying  the  condition 


(II-3) 


For  Rayleigh  scattering  the  four  charactoristi c functions  are 


22 


(u)  = f (1  - (1  + 2u^] 


= fe 


^2(u)  = I (1  - P^) 


'/(u)  = I (1  - p^) 


(II-4) 


For  Rayleigh  scattering  the  scattering  and  transmission  matrices 
are  written  in  terms  of  the  X-  and  Y-functions  as  follows: 


S(p,  (f;  Pq,  (^q)  = Q 


I ■?  (m;  i^o) 


P 1/2  2 V2,.,. 

+ (1  - p'^)  (1  - Pq'^)  S^'Mp,  4>;  Pq,  Oq) 


+ (p,  i>;  pQ,  >Jq) 


(11-5) 
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T iu,  r,  Uq,  ^0^  = Q j T uj 


4 


+ (1  - u)  (1  - u/)  T ^')(p,  o; 


+ T ^ '(p,  o;  Uq,  ({iq) 


where  Q is  as  defined  in  Chapter  I (Equation  1-30),  and,  fo 


T(') . 


(J--  1)  T<*>  = [y,(M)x, (■.:„)  - x,(„)y,(p„)]  T (')(-_ 


A • 
> >•  t 


The  V are  the  azimuth  dependent  components  for  i = 1 , 2 of 
the  Rayleigh  scattering  phase  matrix  given  in  Chapter  I (Equation  1-32 
and  Equation  1-33).  The  azimuth  independent  terms  are 


f(p) 

x(u) 

0 


n/24>(u) 
Vf  C(p) 
0 


t(Po) 


X(Pq) 
V2  c(pq) 


\ 


/ 


/ 


C(u) 
c (v) 
0 


Vz  n (p) 
Vz  e(p) 
0 


0 

0 

0 


Upq) 


0 (pq) 


V2’h(Po)  vre(p.) 


and 


(II“9) 


(?T  - t)  ^ ^®^(p;pq) 


'Pq  VI 
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C(u)  \/2r  (u)  0 


0 


'I'(Mq)  x(mq) 


c (u) 

V2  e(p) 

0 

V2  4)  (Uq) 

V2  4(Pq) 

0 

0 

\ 

0 

0 

0 

/ 

f'{u) 

4>(u) 

'<^o> 

c(Po) 

0 

x(m) 

0 

v'2  n (pq) 

9('Uq) 

0 

0 

0 

V 

(II-IO) 

where  the  T,  n,  a,  and  O are  combinations  of  the 

X3-,  X^-,  Y3-,  and  Y^-functions  and  their  moments.  One  of  the 
differences  between  Chandrasekhar ' s analysis  and  Mullikin's  is  in  the 
exact  form  of  the  combination  of  these  functions.  Chandrasektiai' ' s 
form  of  these  expressions  can  be  found  in  Chandrasekhar  (1950)  and  in 
Chandrusekhar  and  filbert  (1954).  Mullikin's  will  be  given  later  in 
Equation  (11-40). 

Given  these  expressions  for  and  in  terms  of  X-  and  Y- 
functions.  Equation  (11-2)  for  the  X-  and  Y-functions  will  thus  consti- 
tute the  desired  solution  of  the  radiative  transfer  problem. 

I Chandrasekhar  (1950)  developed  an  iterative  metliod  for  their  solution, 

i the  solution  for  small  o[)tical  thickness  as  a first  approximation. 

L 


2f. 


A problem  arises  for  the  case  v/hen 


I 


lij(y)  du  = 1/2 


(II-ll) 


When  this  equality  holds,  Chandrasekhar  recognized  that  solutions  of 
Equations  (II-l  and  II-2)  for  the  X-  and  Y-  functions  are  no  longer 

• (3) 

unique.  For  Rayleigh  scattering  this  condition  occurs  for  'ij' 
(Chandrasekhar's  ).  For  this  "conservati ve"  case  the  solution  is 
given  by  a two-parameter  family.  Chandrasekhar  defined  his  standard 
solutions  by  including  the  additional  constraints 


X(u)'I'(u)  du 


= 0 


Y(y)^^  (u)  dvi 


0 


0 


(H-12) 


which  he  determined  by  a consideration  of  the  flux  equations  at  the  | 

boundaries  of  the  atmosphere. 
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Since  Chandrasekhar  (1950)  first  developed  his  solution  to  this 
problem,  there  have  been  several  publications  showing  different  aspects 
of  the  solution  for  sma 1 1 and  moderate  optical  thickness  (i  ^ 1)  . 
Tables  of  the  X-  and  Y-  functions  have  been  calculated  by  Sekera  and 
Blanch  (1952),  and  Sekera  and  Ashburn  (1953).  Chandrasekhar  and  Elbert 
(1954)  and  Coulson  et  al . (1960),  carried  the  solution  one 
step  further  and  published  tables  of  the  Stokes  parameters  for  radia- 
tion emerging  from  the  bottom  and  top  of  the  atmosphere.  Sekera  (1957) 
has  illustrated  graphically  many  of  the  features  of  the  polarization  of 
the  downward  radiation.  Coulson  (1959)  has  shown  graphs  of  the  upward 
intensity  and  polarization. 

B.  THE  SOLUTION  FOR  LARGE  OPTICAL  THICKNESS  AND  UNIQUENESS 

CONSIDERATIONS 

Chandrasekhar  attempted  to  extend  his  calculation  of  the  solu- 
tion to  larger  optical  thickness  (x  > 1 ) but  was  restricted  by  com- 
putational instabilities.  A similar  attempt  by  Blanch  was  restricted 
by  computer  size.  This  type  of  extension  was  later  attempted  by 
Sekera  using  the  solution  for  t ->  «■  as  a first  approximation  rather 
than  that  for  small  x . However,  some  of  the  functions  failed  to 
converge  sati sfactori ly  even  after  many  iterations  and  it  was  noted 
that  the  iterated  values  were  of  an  oscillatory  nature.  Thus  it  was 
apparent  that  modifications  had  to  be  made  to  Chandrasekhar's  method 
of  solution  for  the  extension  of  results  for  larger  optical  thickness. 
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Mullikin  (1962a,  b,  c,  d;  1964a,  b;  1965)  undertook  an  exten- 
sive study  of  the  uniqueness  problem.  He  determined  that  the  diffi- 
culty arose  in  the  usual  heuristic  reduction  from  the  transfer  equa- 
tions to  the  equations  for  X-  and  Y-  functions.  This  process  ignored 

some  equations  and  this  loss  of  information  introduced  extraneous  solu- 

(3l 

tions.  This  had  been  recognized  in  the  case  of  by 

Chandrasekhar,  but  Mullikin  realized  it  also  extended  to  and 

(21 

\p'  Mullikin  developed  additional  equations  in  the  form  of  linear 
constraints  which,  along  with  Chandrasekhar's  X-  and  Y-  equations,  do 
select  unique  S and  T matrices.  He  determined  these,  as  discussed 
below,  by  relating  the  integral  equations  back  to  the  original  radiative 
transfer  equation.  These  constraints  contain  information  that  is  lost 
in  going  from  the  i ntegro-di fferential  equations  to  the  coupled  non- 
linear equations . 

Mullikin  first  considered  the  work  of  Busbridge  (1960)  who 
demonstrated  the  existence  of  solutions  to  the  X-  and  Y-  equations 
(Equations  (II-l)  and  (II-2)).  Busbridge  investigated  the  auxiliary 
integral  equation 


1 


J(x.p)  = exp  (-J-)  + j exp  (-^^^^)J(y,p)dyd\ 


(11-13) 


where,  for  isotropic  scattering,  J is 
showed  that  a solution  to  Equations  Il-l 


the  source  function.  She 
and  II-2  is  given  by 


29 


I 


X„(u)  • J(0,u) 
Yq(m)  = J(t,u) 


(11-14) 


and  that  and  are  defined  for  all  complex  p,|u|>0  , to  be 
real  and  non-negative  for  u ^ 0 and  to  be  analytic  in  the  extended 
complex  ..-plane  except  near  u = 0. 

Mullikin  (1962a,  b,  c,  d),  following  her  proof  of  existence  of 
solutions,  turned  to  the  question  of  uniqueness,  since  he  realized 
that  there  appeared  to  be  a multiplicity  of  solutions.  He  returned 
to  the  transport  equations  (e.g..  Equation  1-27)  and  noted  that  for 
the  physically  significant  solution  to  most  scattering  problems,  the 
X-  and  Y-functions  must  have  extensions  to  complex  u as  analytic 
functions  in  'ur’O  . Following  Busbridge,  he  considered  the  exten- 
sion of  Equations  II-l  and  II-2  to  the  complex  plane,  by  rewriting 
them  in  the  form 


1 


X(v)c (v) 

v+3- 


0 


* n})^/  = 1 

0 


(11-15) 
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with  ^ any  complex  number  not  in  the  interval  [-1,  1].  As  a sys- 
tem of  linear  equations  in  X(j.)  and  Y(^),  these  equations  have  a 
unique  solution  where  the  determinant  does  not  vanish.  This  determin- 
ant is  the  function 

' 0 

(11-17) 
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Equations  (11-15)  and  (11-16)  can  then  be  rewritten 


1 1 


(11-20) 


Mull i kin  has  shown  that  these  equations  define  a meromorphi c extension 
of  X and  Y to  the  complex  domain  lj-|>0  , giving  functions  analy- 
tic in  |^|>0  except  for  possible  poles  at  the  zeros  of  A . 

Following  Busbridge  (1960),  Mullikin  considered  the  character- 
istic equation 


(11-21) 


The  uniqueness  or  multiplicity  of  solutions  to  the  X and  Y equa- 
tions will  depend  upon  the  roots  of  this  equation,  which  will  in  turn 
depend  upon  the  scattering  characteristic  functions  i]j(v),  given 
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for  Rayleigh  scattering  by  Equation  (II-4).  Mullikin  shovjed  that  four 


cases  can  occur 

I \ has  no  zeros 

II  The  only  zeros  of  X are  at  1/k  v;here 

0 < k < 1 (the  non-conservative  case) 

III  The  only  zeros  of  X are  at  + 1/k  where 

k = 0 (the  conservative  case) 

IV  X has  a zero  at  + 1 . 

Mullikin  assumed,  for  simplicity,  that  case  IV  never  occurs,  which  is 
true  for  Rayleigh  scattering. 

For  case  I,  X with  no  zeros,  Mullikin  (1962a)  showed  that 

the  only  bounded  solutions  to  Equations  (11-19)  and  (11-20)  (and 

hence  to  Equations  (II-l)  and  (II-2))  are  given  by  the  and 

solutions  of  Busbridge  (Equation  (11-14)).  This  is  the  case  for 
(4) 

of  the  Rayleigh  scattering  problem  (see  Equation  (1 1-4)). 

For  case  II,  the  non-conservative  case  with  zeros  of 
X at  1/k  , Mullikin  (1962a,  d)  showed  that  all  solutions  to 
Equations  (Il-l)  and  (II-2)  are  given  by 

(11-22) 
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where  the  a and  b are  constants  constrained  only  by  the  condition 
(a^  - b^)  (a^  - B^)  - 2a  ak  - 2Bbk  = 0,  and  where  and  are 

the  solutions  of  Busbridge  (Equation  11-14).  Thus,  vvc  liave  a family 
of  solutions  depending  upon  the  values  of  the  parameters  a and  b. 
Mullikin  (1962a,  d)  considered  the  behavior  of  these  solutions  at  the 
zeros  of  A to  derive  the  additional  constraints. 
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When  there  is  such  a multiplicity  of  solutions,  then  one  must  deter- 
mine which  solution  corresponds  to  the  physical  problem.  For  most 
scattering  problems  the  solution  that  is  analytic  for  all  > 0 is 
the  desired  solution.  This  analytic  solution  can  be  selected  from  all 
the  possible  solutions  by  demanding  that  the  system  of  Equations  (11-19) 
and  (11-20)  be  well  behaved  at  the  zeros  of  X . Thus  when  X *■  0 
tlie  R.fl.S.  of  Equation  (11-19)  and  (11-20)  must  also  go  to  zero.  This 
imposes  the  constraints  that 
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and 


/ 


l-kv 


Y(v)'>^(v)J  _ 


(11-30) 


, f X(v)iJ;(v) 

' J 1+kv 


1 


-kv 


0 


0 


(11-31) 


By  comparison  with  Equation  (11-27)  and  (11-28)  we  see  that  this  imposes 
the  constraint  that  a = b = 0.  Then,  by  substituting  a = b = 0 into 
Equations  (11-22)  and  (11-23),  we  see  that  the  solutions  corresponding 
to  the  physical  problem  in  this  case  are  the  analytic  solutions, 
and  Yq  of  Busbridge.  Case  II  is  encountered  in  the  Rayleigh- 
scattering  problem  for  and  (Equation  (II-4)).  This  non- 

uniqueness was  not  recognized  previously. 

Finally  we  consider  Case  III,  the  conservative  case  (iji^^^) 
where  k = 0 (the  zeros  of  X are  at  «).  Here  Mullikin  (1962d) 
showed  that  all  solutions  to  Equations  (II-l)  and  (II-2)  are  given  by 
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x(u)  = X^(m)  + au(X^(u)  + Yq(u)) 

+ bpLlXpCu)  + m(Xq(m)  + VQ(y))J 


where 


Y(vj)  = Yq(u)  - ay(Xp(u)  + Yq(u)) 

- by[YY  (u)  - u(Xo(p)  + Y^(u))] 


Xl  + Xl 

—JT 


X = / V ij;(v)Xjj(v)dv 


= J"  v%(v)Yq(v)c1v 


with  n = 0,1  and  the  a and  b are  constants  constrained  now  by 

b(by^  2ay  - 2)  = 0. 


i 

J 


1 1 
2, 


2a  / v^4)(v)dv  = X ^ Y ( j )■;■( v)dv 
0 0 


(11-37) 


y*[X(v)  - Y(v)]v';  (v)d', 
0 


(11-38) 


Once  again  the  analytic  solution  will  have  a = b - 0,  with 
the  solutions  to  Equations  (11-19)  and  (11-20)  being  given  by  and 
Yq  . In  this  case,  however,  for  Rayleigh  scattering,  Mullikin  (1966) 
has  shown  that  the  analytic  solution  is  not  ttie  desired  solution.  The 
azimuth-independent  solution  is  not  the  desired  solution.  The  azimuth- 
independent  parts  of  the  Rayleigh-scattering  transmission  and  scatter- 
ing matrices  are  deteniiined  from  the  X3  and  Y3  solutions.  There- 
fore, Mullikin  examined  the  azimuth-independent  part  of  the  source 
function  given  by 
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1 1 

J(x.u,^  ""  J J exp  J(y.v,^)dy  ^ 

0 0 

+ p(u.>)  exp 

(11-39) 


This  does  not  go  to  zero  as  >.  -*■  ® but  depends  quadrat! cal ly  on  ^ . 
Thus,  in  Equations  (11-32)  and  (11-33)  the  a and  b must  be  non-zero. 
Mullikin  (1964)  derived  linear  constraints  based  on  the  original  trans- 
fer equation  which  ensures  the  required  quadratic  behavior  of  the  solu- 
tions as  3- 

Sekera  (1966a,  b)  showed  how  to  express  the  solution  of  the 
Rayleigh-scattering  radiative  transfer  problem  in  terms  of  combinations 
of  the  analytic  X's  and  Y's.  The  linear  constraints  derived  by 
Mullikin  are  included  in  the  constants  used  in  the  combining  of  the 
analytic  X's  and  Y's  to  determine  the  azimuth-independent  part  cf 
the  scattering  and  transmission  functions.  Following  Mullikin,  Sekera 
showed  the  following  relationships  bctvjeen  the  K-  and  L-  functions 
(analogous  to  Chandrasekhar ' s ik  . ^ , etc.,  functions  of 

Equation  (11-9)  used  for  the  azimuth-independent  part  of  the  solution) 
and  the  X^-  , X^-  , Y^-  and  Y^-  functions. 
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= 'VM  = (c  + by)u  x^{u)  + (a  + bp)u  Y3(;:) 
L^(u)  = ^(u)  = (-  a + by)u  X3(u)  + (-  c + by)y  Y^Cy) 


K,(y)  = (D(y) 


= (1  + c'y  + b'y^)  X3(y)  + (a'  + b‘y)y  Y3(y) 


L2(u)  =n(y)  = (-  a'  + b'y)y  X3(y)  - (1  - c'y  + b’y'^)  Y3(y) 

l<3(y)  = xM  = (1  + ey  - fy^)  X^(y)  + (-  g + fy)y  Y^(y) 

l-3(y)  (y)  = (g  + fy)y  X^(y)  + (1  - ey  - fy^)  Y^(y) 

2K4(y)  H 2;(y)  = (c  + f y)y  + (a  - f y)y  Y^(y) 

2L^(y)  H 20(y)  = (a  + f'y)y  X4(y)  + (-  c + f’ y)y  Y^ty) 


whe  re 


b = - ^ A 


b'  = 2mb 


f’  = ]■  A' 


f = 2mf 


= - br  + Gf  a'  = - b'r  + 2(-  s + nG)f 

= br  + Gf  c'  = b'r  + 2(-  s + nG)f 


e=-2(r+nr)b-Gf 


g = - 2(r  + nr)  b + Gf  (11-43) 


2A 


1 

m + sr 


2A' 


1 

m - rG 


P _ ■ “3 

2(B]  - 6‘)  - T(e'  + 6^) 

^0"^2 

2(a'  + Qt^ 

(11-44) 

r = a^aj  - 

" “o^o  " 

s = 6q6]  - 6^6q 

“O  " ^ ■ T 

^0  1 1 "’O  ^^^3^ 

0*1  = f mi  [P3] 

3i  = f mi  [Q3]  (i  = 1.  2 .••) 

“0  ~ ^ " 8 "’0  ^^4^ 

^0=1"! 

“i  = l'"i 

= I"’!  [Q4]  (i  = 1.  2 •••) 

P,  = Xi(y)  + Yi(y) 

Ql  = Xi(p)  - Yi(u)  (i  = 3.  4) 

(11-45) 
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and  m^Lf]  stands  for  the  i-th  moment  of  the  function  F,  i.e.. 


4 


1 

[F]  = J F(x)  x’  dx  (11-46) 

0 

These  relationships  only  hold  for  the  case  wq  = 1 , for  perfect 
Rayleigh  scattering  with  no  absorption.  For  Wq  > 1 , the  K-  and  L- 

functions  cannot  be  expressed  by  the  simple  relationships  above 
(Equation  (11-40)).  The  reduction  of  the  singular  integral  equations 
for  the  K-  and  L-functions  is  then  much  more  complicated,  and  will  not 
be  discussed  here. 

This  is  the  case  where  Chandrasekhar  (1960)  recognized  the  non- 
uniqueness and  specified  his  "standard"  solutions  by  consideration  of 
the  flux  equations.  If  one  would  like  to  use  Chandrasekhar's  solution 
it  is  necessary  to  introduce  the  moment  conditions  for  the  "standard" 
solutions  in  the  equations  above.  Since  the  "standard"  solution  is 
defined  by  the  relations 


(11-47) 


one  has  to  substitute  in  Equation  (11-44)  1/:  0 in  oi  d • .'I  t.iin 

the  corresponding  forms  of  the  equations  fo>  "stand  h r • 

This  leads  to  A = 0,  h 0,  b'  0 , but  w*’  scf  tbi' 


•i 

1 


1 
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); 


Equations  (II-5)  through  (II-IO)  along  vnth  (11-40)  thus  con- 
stitute the  complete,  unique  solution  to  the  Rayleigh-scatteri ng  radia- 
tive transfer  problem  for  all  optical  thicknesses. 

Mullikin  (1962d)  also  obtained  from  the  constraints  he  developed 
a new  formulation  of  the  equations  that  was  more  suitable  for  numerical 
computations.  These  were  particularly  useful  for  atmospheres  of  large 
optical  thickness.  His  computational  scheme  will  be  outlined  in  the 
following  section. 

Following  these  theoretical  considerations  by  Mullikin  and 
Sekera,  a computer  program  was  developed  to  compute  the  X-  and  Y-func- 
tions  and  the  K-  and  L-functions.  Carlstedt  and  Mullikin  (1966)  calcu- 
lated and  published  tables  of  X-  and  Y-functions  for  the  case  of  iso- 
tropic scattering,  which  has  the  single  characteristic  function  ^ 
and  only  one  set  of  X-  and  Y-functions. 

Sekera  and  Kahle  (1966)  extended  the  computation  to  the  case  of 
Rayleigh-scattering.  They  published  tables  of  the  X-  and  Y-functions 
and  the  K-  and  L-functions  for  a complete  range  of  values  of  i and 
n . These  tables  are  reproduced  here  in  the  Appendix.  Their  calcula- 
tion and  application  arc  the  subject  of  the  remainder  of  this 
dissertation. 
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C.  METHOD  OF  COMPUTING  THE  X-  AND  Y-FUNCTIONS 

Mullikin  (1962c,  d)  and  Carlstedt  and  Mullikin  (1966)  have 
shown  how  to  solve  for  the  X-  and  Y-functions  for  the  three  cases  dis- 
cussed in  the  previous  section.  The  portion  of  their  discussion  rele- 
vant to  the  computer  calculation  of  the  X-  and  Y-functions  for  Rayleigh 
scattering  is  outlined  below,  in  sufficient  detail  to  enable  the  reader 
to  develop  computer  programs  to  perform  such  calculations. 

CASE  1 --  Unique  Case 

We  will  look  first  at  the  unique  case,  \ with  no  zeros.  We 
define  the  following  functions: 


e(t)  = I tan'^ 


.t)l 

. ^0 

ttjj 

(11-49) 


N(z)  = exp 


/ 


e(t) 

t - z 


dt 


(11-50) 


where 


= 1 


2^.2  J ^ In  (11-51) 


i 
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and  final  ly 


N(0)  = 


(11-52) 


Then  Chandrasekhar's  H- function  (Chandrasekhar,  1960)  can  be 
shown  to  be  related  to  the  N-function  (Mullikin,  1964a) 


r ^ 

-1/2 

r ^ 1 

(0)'  = 

1 - 2 J"  ip{t)dt 

exp 

. /'e(0 

/ t+z 

0 

*/ 

^ 0 

(11-53) 


Mullikin  then  shows  that  the  X-  and  Y-functions  can  be  written 
for  Y outside  [-1,0] 


X(z)  = H(2) 


1 


1 + 


(11-54) 


Y(z)  = H(z) 


g(z)  + z 


t)[q 

i ^ J 

- g(z)] 

Hi 

JiTaI 

dt 


(11-55) 
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where 


A(t)  = ^ 


(11-56) 


and 


f = £-^  g 


= P-  q 
2 


(11-57) 


p(z)  = - L(p)(z)  + 


(11-58) 


q(z)  = L(q)(z)  . 


(11-59) 


and  the  operator  L is  defined  as 


L(p)(z)  = 


z exp  (-  t/z) 
■ H(z) 


1 

/ 


>P  (t)p(t)  dt 

H(t)  \Mt)  (t  + zT 


(11-60) 


The  equations  for  p and  q are  Fredholm  equations  in  the  interval  of 
interest,  and  can  be  solved  by  iteration.  The  convergence  of  these 
iterations  is  very  rapid,  and  is  faster,  the  larger  the  value  of  t , 
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The  rest  of  the  problem  is,  then,  just  quadratures.  This  method  is 
good  for  all  values  of  t . 

CASE  2 --  rion-unique.  Non-conservative  Case 

For  this  case  the  X-  and  Y-functions  can  be  written  as 


X(z)  = H(z)  = 


1 + 


kzN(O)  f(z)  - f(l/k) 
1 - kz 


(11-61) 


+ z 


f HUlf.m  - .fc.zjl.. 

J H(t)  A(t)  (t  - z) 


dt 


Y(z)  = H(z) 


(11-62) 


where  f and  g are  defined  in  terms  of  p and  q as  before,  but  p 
and  q are  more  complicated. 


p = h^  - ch2  q = h^  - dh 


3 '^"4 


(11-63) 
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where 


(-  1) 

= - L(h2)  + zN(-  z)  exp  (-  |) 

*^3  ^ N(o)^ 

^4  = L(h4)  + zN(-  z)  exp  (-  j) 


(11-64) 


and  the  operator  L is  now 


These  are,  again,  rapidly  converging  Fredholm  equations. 
The  constants  c and  d are 


h^(l/k)  + h^(-  1/k)  exp  (-  ki) 
h2(l/k)  - h2(-  1/k)  exp  (-  k j 


(11-66) 


h3(l/k)  - h^(-  1/k)  exp  (-  ki) 
^ ” h^(l/k)“-  h^(-  1/k)  exp  (-  ktT 
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These  results  are  derived  from  allowing  k to  approach  zero  in 
the  above  case.  Now 


X(z)  = H(z)  1 + [1  - zN(0)]  f(z) 


(11-67) 


Y(z)  = H(z)  [1  - zN(0)]  g(z) 


N(0)d  ^ ^ i 

1 

[9(t] 

- g(z)]  ! 

-y-  Z + Z y 

H 

Ct  - zT^M 

where  now  the  constant 


2 

1 .M3 

nT^ 

1 

T + 2 

1 - / e(t)  dt 

2 M 

nW  ”4 

0 

(11-69) 


and  the  constant  c = 0 , and 


/ 


ip(t)(l  - kt)h-(t) 
H(t)  A(t) 


dt 


(11-70) 


Thus  all  the  values  of  the  X-  and  Y-functions  have  been  re- 
duced to  the  solution  of  rapidly  converging  Fredholm  equations  and 
quadratures.  Once  the  X^-,  X^- , Y^-,  and  Y^-functions  have  been 
calculated,  the  K-  and  L-functions  are  quickly  calculated  using  Equa- 
tions 11-40  through  11-46  in  Section  B of  Chapter  II.  From  these 
functions,  the  solution  to  several  problems  in  radiative  transfer  can 
be  quickly  determined.  A description  of  hovj  these  functions  can  be 
used  is  given  in  the  next  chapter. 
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III. 


DESCRIPTION  OF  THE  TABLES  AND  THEIR  USE 


In  this  chapter  the  tables  of  X-  and  Y-  and  K-  and  L-functions 
which  are  given  in  the  Appendix  are  described.  The  method  of  finding 
the  Stokes  parameters  from  these  functions  is  written  out  in  detail. 
Other  uses  for  the  tables  are  also  suggested.  Finally,  the  accuracy  of 
the  tables  is  discussed. 

A . DESCRIPTION  OF  THE  TABLES 

Table  1 in  the  Appendix  contains  the  X-  and  Y-functions  for  the 
: values  0.15,  0.25,  0.50,  0.70,  1.0,  2.0,  4.0,  8.0,  16.0,  and  100.0. 

The  functions  are  listed  for  all  „ values  from  0.00  to  1.00  in  steps 
of  0.02.  Table  2 contains  the  moments  of  the  X-  and  Y-functions  of  or- 
der zero  througli  four,  for  the  same  ^ values  as  in  Table  1. 

Tables  3 and  4 give  the  K-  and  L-functions  and  their  moments  for 
the  same  parameters.  All  values  are  given  to  five  figures.  The  deci- 
mal part  of  the  number  is  followed  by  an  E and  the  power  of  10  to  wliich 
the  decimal  part  should  be  raised;  i.e.,  0.17452  E-01  = 0.017452. 

For  the  largest  value  of  the  optical  thicknc.ss  given  here, 
t = 100,  all  the  X-  and  Y-functions  except  and  Y^  have  essentially 
reached  their  values  for  an  atmosphere  of  infinite  thickness  (the  X- 
functions  approaci)  the  H-f unctions  [Chandrasekhar,  1950]  and  the  Y- 
functions  approach  zero).  In  the  conservative  case  (X^  and  Y-)  the 
values  of  the  X-  and  Y-functions  approach  their  limiting  values  as  1/:  , 
rather  than  exponentially  as  in  the  other  cases.  Since  the  values  of 
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the  K-  and  L-functions  depend  upon  and  Y^,  they  too  have  not  yet 
reached  their  value  for  an  infinitely  thick  atmosphere. 

B.  CALCULATION  OF  THE  STOKES  PARAMETERS 

For  a plane-parallel  atmosphere,  Rayleigh-scattering  atmosphere, 
— ► 

the  Stokes  parameters  I = (I^,  I^,  U,  V)  of  the  radiation  emerging 
from  either  the  top  or  the  bottom  can  be  written  as  the  sum  of  three 
terms : 


-*•  ■'■f  nl  '"n  1 

I(t;  p,({))  = I'  ^(x;  p,  Pq)  + r ^(i;  p,  Pq,  (p  - (p^) 


+ p)rT:, 


(ni-1) 


r ^[t;  p^  Pq,2(!?)  - (Jq)] 


where,  in  the  same  manner  as  with  the  scattering  and  transmission  m,a- 
triccs  and  the  Rayleigh-scattering  pliase  function,  the  first  term  is 
azimuth-independent,  and  the  second  and  third  terms  contain  the  elements 
with  tlie  cosine  or  sine  of  (o  - f>Q)  or  2(p  - ijp)  , respectively.  For 
Rayleigli  scattering  is  zero,  as  are  all  three  V terms,  so  that 

►(i ) 

I'  can  be  reduced  to  a two-  or  three-element  matrix.  In  Equation 
(111-1),  (n,  i>)  and  (..q,  .-q)  are  the  sets  of  the  usual  directional 
parameters  of  the  emerging  radiation  (.. , ;■)  and  of  the  external 
parallel  solar  radiation  (,.q,  ;q)  irradiating  the  top  of  the  atm.osphere. 
For  the  upward  radiation,  x = 0 and  is  positive;  for  the  downward 


1 

i 
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radiation,  r = (the  total  optical  thickness  of  the  atmosphere) 

and  M is  negative.  If  -nf  denotes  the  net  flux  of  the  external  ra- 
diation,  then  by  combining  Equations  (1-39)  and  (1-40) (defining  I in 
terms  of  S and  T ) with  Equations  (11-5)  through  (II-IO)  and  (11-40) 

^ — ► -i — ► 

(defining  S and  T in  terms  of  the  X-,  Y-,  K- , and  L-functi ons ) , the 
azimuth-independent  tenns  can  be  written  in  the  form  [Chandrasekhar, 
1950;  Sekera,  1963,  1966b]: 


+ p.  - Un)  = 


M,  Uq) 


lj°)(0;  u,  Uq) 


3_  ^0^0  r7,  , . 

16  Uq  + u ^ 


K (Pn)  - L (p)  • L (Pfi)  • I" 


- u,  - Pq)  = 


U,  Pq) 


P.  Pq) 


(III-2) 


•3  Pn‘*^n 

16  Uq  - p ‘ ^ ^‘'0^ 


L (p)  • K (Pq)]  • E 


(III-3) 
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I 


i 

! 

I 

v;here  the  two-by-two  matrices  k and  L have  the  fonii  i 

N^(u)  , 

H^{u)  , V^N^(u) 

The  functions  are  given  in  Table  3 as  functions  of  u for  a 

particular  value  of  the  parameter  and  for  u)g  = 1 . 

Assuming  the  unpolarized  external  irradiation  = 1/2Fq, 

then  the  azimuth-dependent  terms  are,  for  upward  radiation, 


(III-4) 

(ri  - K,  L) 


(III-6) 


where  C = (3/32)lvq;.qFq 

Pq)  = X.(-m)  X.(wq)  - Y.(u)Yi(Po)  (i  - 1,2)  (111-7) 
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For  downward  radiation, 


-*■(] ) 

I (t-|  , - y , - Pq,  qj  - (|ip)  = 


4Cm 


0 ? ''/2  _ 1/2 

“ 0 - (1  - u^) 


Mq  - U 


I ' “^1  ’ ■ ■ 1^0’  ^ ” ^0^ 

(1) 


Ip  (T]  ; - M,  - Mq,  $ - $p) 

/ 


.. 


u cos  (p  - $q) 


u,  Uq) 


\ 


- sin  (0  - (p  ) 


(III-8) 


/ \ 


•>  ( 2 ) / 

j (t^  ; - M , - Uq,  2(<()  - 4>-) ) = 


\ 


(2) 

Ip  (fp  - ^Jl  - tJQ,  2{<p>  - iJ)q)) 
- p,  - M„,  2(^  - 


/ 


~ 0 - uh 


I - cos  2(t  - ,>jj)  \ 


COS  2(4>  - ((.q) 
2 n Sin  2(4>  - 4>q) 


; M.  pq) 


/ 


(III-9) 


I 

i 

I 

, 

I 


I 
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p 


where 


p,  Uq)  = X.(u)  Y.(uq) 

- Y.(m)  X.(uq)  (i  = K 2) 


(IIMO) 


Illustrations  of  the  use  of  the  tables  to  find  the  intensity  of  radia- 
tion emerging  from  the  top  and  bottom  of  a plane-parallel  Rayleigh- 
scattering  atmosphere  will  be  given  in  Chapter  IV,  Section  A. 

C . ADDITIONAL  USES  OF  THE  TABLES 

The  primary  purpose  of  Tables  1 through  4,  as  discussed  above, 
is  to  provide  numerical  values  of  the  functions  needed  to  compute  the 
intensity  and  polarization  parameters  of  the  radiation  emerging  from 
the  top  and  from  the  bottom  of  a plane-parallel  planetary  atmosphere 
with  Rayleigh  scattering.  Equations  (III-l)  through  (III-IO)  give  the 
computational  scheme  for  this  application. 

In  addition,  the  tables  can  be  used  to  compute  the  following 
quantities,  mentioned  in  the  five  headings  below. 

Chandrasekhar's  functions  Y£(’.  u).  m)«  s(t) 

These  are  needed  to  compute  the  effect  of  ground  re- 
flections governed  by  Lambert's  law  on  the  average  intensity  and  net 
fluxes  of  the  emerging  radiation.  Lambert  reflection  is  where  the  sur- 
face reflects  unpolarized  light  uniformly  in  all  directions,  independent 
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of  the  direction  and  polarization  of  the  incident  light.  As  shewn  by 
Sekera  (1966b),  these  functions  can  be  computed  from  the  values  of 
the  functions  , L^-  and  their  zero  and  first  moments.  Using  the  mo- 

ment notation  of  Equation  (11-46),  we  can  write  [Sekera,  1966b,  Eq. 
(150),  p.  51] 


Y(>(y)  = i!-i[KT(T,  u)  + L,(t,  u)]  + 2?.j[K,(t,u)  + L,(t,  u)] 


(III-ll) 


Y^(y)  = 


£^[K3(t,  p)  + 13(1,  p)]  + + L^v':,  ;.)] 


(111-12) 


+ K3  + L.|  + L3]  - 2.e2ni^[K2  -f 


K4  ■'■1-2  + L^) 


whore 


(111-13) 


l^i  = g vO  + L.^2(i  , p)]  (i  - 1 , 2)  (111-14) 


If  ground  reflection  is  according  to  Lambert's  law  with  reflec- 


tivity A,  the  terms  that  must  be  added  to  the  intensity  vectors 
+ M.  - Uq)  and  - u,  - Uq)  have  the  form 

[Chandrasekhar,  1950  ,p.  279),  for  upward  radiation 


I (0;  u,  - Mq)  = 


4[1  - A?(t)] 


Uq)  + (t. 


Yj^(t,  u) 


Yp(T,  p) 


and  for  downward  radiation. 


^pQ^O 

I (Yii  - p.  - Pn)  = r [yci(^  > Pn)  ] 

' 0 U 0 r 0 J 


1 - Y £ ( t > P ) 


1 - Y^(t,  u) 


(III-16) 
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This  is  defined  as 


1 2tt 


J(Oi  Uq)  - 


0 0 


(III-17) 


or 


1 2tt 


h f f - M,  4))  + I^(t^;  - u,  ^))]  d-4  di 


0 0 


(111-18) 


Upon  integration  with  respect  to  4>  , these  reduce  to 


ln,„[l<“>(0;  * u,  - i,„)  t - P„)]  (UM9) 
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J 


or 


[ 


1 

1 


- p. 


Mq) 


Po>] 


(III-20) 


respectively.  Sekera  has  shown  [1966b,  Eq.  (14Q),  p.  47]  that  these 
expressions  can  easily  be  computed  in  terms  of  the  K-  and  L-functions 


4 *• 

J(0,  Uq)  = (V4)[^K^(t^,  Pq)  - 2]  (III-21) 

n=l 


and 


4 

J(T,.  Pq)  ■ 

n=1 


>)  - 2e 


(III-22) 


(c)  Global  Radiation,  and 

The  downward  global  radiation,  G^  , is  the  total  ra- 
diation reaching  a unit  surface  area  of  a planet,  and  includes  both  the 
diffuse  sky  radiation,  H^,  and  the  direct  (but  attenuated)  solar  ra- 


I 


Deirmendjian  and  Sekera  (1954)  give  the  following  expressions 
for  downward  radiation: 


S = 


(HI-23) 


i 

i 

i 

I 


^d  "" 


Yo(hn)  + 


P vhq;  ■ . r^‘-Q^  „ ^0 


-T/Pf 


2(1  - As) 


(III-24) 


- Hd  + 5 tiFqI.q 


Yp(Uo) 


0' 


2(1  - As) 


(111^6) 


The  upward  radiation,  as  given  by  Coulson  (1959),  is 


1 - 


Y„(uo)  ^ 1- 


1 - As 


(HI-26) 


6? 


of  reflection 


The  intensity  matrices  of  the  contribution  to  the 
emerging  radiation  from  ground  reflection  governed  by  a more  general 
lav/  than  Lambert's  (see  (a)  above)  can  be  expressed  in  a rather  com- 
plicated form  containing  the  reflection  and  transmission  matrices. 

Hence,  after  a lengthy  reduction,  we  obtain  the  expressions  for  these 
intensities  that  contain  the  functions  (i  - 1,  2,  3,  4),  as 

well  as  X.,  Y.  (j  = 1,  2).  These  expressions  for  Fresnel's  law 

J sJ 

(specular  reflection  with  partial  linear  polarization)  can,  for  example, 
be  found  in  papers  by  Sekera  and  Frazer  (1953)  and  Frazer  (1965). 

( e ) Computation  of  the  characteristics  of  the  internal 

radiation  field 

The  intensity  and  polarization  of  the  diffuse  radiation 
at  any  level  within  the  atmosphere  can  be  computed  from  the  principles 
of  invariance  that  express  the  intensity  vector  of  the  upward  and  down- 
ward radiation  in  terms  of  the  diffuse  reflection  and  transmission  by 
the  layers  above  and  below  the  reference  level.  Using  the  diffuse  re- 

•4— ► 4— ► 

flection  and  transmission  matrixes  S and  T defined  in  Equations 
(1-39)  and  (1-40)  respectively,  we  have  for  the  upward  radiation  at 
the  level  i (0  < i < t-|  ) 
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uI(t;  + m,  4>)  = 4 S (ti 


-»■  -t/Uq 

U»  ^Q’ 


(S(t^ 


t;  u,  $;  m'  , ^)' ) • uI(t;  , ({)')} 


(111-27) 


and  for  the  downward  radiation 


^ 1 — ► 

Vil(T-,  - p,  4>)  = ;^T(t;  p,  ());  Pq,  (},q)  • F 


(S(t; 


) • pI(t;  + p,  $)} 


(III-28) 


where  the  symbol  used  for  hemispherical  integration  is  defined  as 


2ti  1 

^ X*  dij  * 

{F(p‘ . (J.’)}  = h J J ^ (III-29) 

0 0 
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These  relations  follow  directly  from  the  physical  meaning  of  the  dif- 
fuse reflection  and  transmission  by  the  sublayers  above  and  below  tiie 
reference  level.  These  equations  can  be  expanded  as  before  into  two 
equations  for  the  azimuth-independent  terms,  and  one  for  each  azimuth- 
dependent  term.  These  integral  equations  must  be  solved  by  successive 
iteration;  to  compute  these  iterations  for  the  azimuth-independent 
terms,  we  need  the  K^-  and  L.-functions  for  the  optical  thicknesses 
T and  T.J  - T , and  for  the  azimuth-dependent  terms  the  functions 
X^,  Y.  (i  = 1,  2). 

The  foregoing  list  of  possible  uses  of  tabulated  functions  is 
by  no  means  complete.  The  applications  mentioned  should  be  regarded  as 
selected  examples.  One  can  easily  extend  this  list  to  include  other 
applications.  For  example,  the  tables  can  be  used  in  problems  of  mixed 
scattering  in  which  the  phase  matrix  is  the  sum  of  a Rayleigh  scatter- 
ing matrix  and  an  additional  matrix  representing  contributions  from 
isotropic,  neutral  (molecular  anisotropy,  resonance  scattering),  or 
aerosol  (turbid  atmosphere)  scattering.  Hov/ever,  some  of  these  appli- 
cations cannot  be  carried  out  yet,  as  the  explicit  formulas  required 
have  not  been  developed. 

D,  ACCURACY  OF  THF  TABULATED  VALUES 

There  are  several  relationships  that  should  be  satisfied  by 
the  functions  , Y^  (i  - 1,  2,  3,  4)  or  their  moments,  and 

v;e  have  used  them  to  check  the  accuracy  and  internal  consistency  of  the 
computed  values  of  these  functions. 
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From  the  physical  fact  that  the  neutral  lines  (i.e.,  the  lines 
where  O'”  Q " 0)  must  pass  through  the  zenith  or  the  nadir 

at  45°  to  the  sun's  vertical  (i.e.,  Q for  J - 0^  = 45°  must  approach 
zero  for  p -*•  1 ) , it  follows  that  the  K-  and  L-functions  for  u = 1 
must  satisfy  the  condition 


Ki  - K3  = K2  - K4  = Li  - L3  = L2  - L4  = 0 (111-30) 


It  can  be  seen  that  all  five  significant  figures  of  the  tabu- 
lated values  satisfy  this  condition  for  every  x . This  accuracy 
may  be  compared  with  the  accuracy  of  previous  computations 
[Chandrasekhar  and  Elbert,  1954;  Sekera  and  Blanch,  1952],  as  shown  in 
Table  A. 

The  relations  for  the  moments  of  the  K-  and  L-functions  can  be 
derived  from  the  linear  constraints  on  the  solution  of  the  integral 
equation  for  the  azimuth-independent  tenns  of  the  reflection  and  trans- 
mission matrices  [Sekera,  1966b,  p.  49].  If  we  introduce  the  following 
abbreviations  for  the  moments 


k) 


(n)  _ 3 


Ki+2]' 


(n)  _ 3 


'-i+2'' 


(i  = 1,  ?-) 


(III-31) 
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then  these  relations  have  the  form 


= kj°^  + = 1 

(111-32) 

M2  H 2k^°^  + = 1 

(111-33) 

M3  E 2[k{^^  - - t[1  - k}°^  + = 0 

(111-34) 

M4  ^ 4[k^’^  - t[l  - 2kJ°^  + 2£^°^]  = 0 

(111-35) 

The  left  side  of  Equations  (III-32)  through  (III-35)  are  tabulated  for 
different  values  iti  Table  B.  The  larger  deviations  from  zero  for 
large  values  of  i result  from  the  loss  of  significant  figures  in  the 
checking  process  rather  than  from  the  inaccuracy  of  the  mom.ents.  The 
moments  of  the  X-  and  Y-functions  must  satisfy  two  sets  of  relations. 

For  the  characteristic  functions  of  the  form 

= a + bu^  (111-36) 

we  have  the  relation  [Chandrasekhar,  1960,  p.  189] 
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I 

I 

1 

-3 


J 

i 


(111-37) 


where  a , g stand  for  the  n-th  moments  of  the  functions  X and  Y , 
n n 

respectively.  Since  a = - b = 3/4,  for  Y^  and  a = - b = 3/8 
for  X^,  Y^  , we  have 

N,  H (»3)„  - I [(c3)2  - (Sj)^  - («3)2  * ■ 1 (IU-38) 

"2  - (6,)^  * (6,)?]  • 1 (UI-39 

where 


stand  for  the  n-th  moment  of  the  functions  X^  and  Y^  , respectively. 
The  left  sides  of  Equations  (III-38)  and  (III-39)  are  tabulated  in 
Table  C. 

Another  set  of  relationships  follows  from  the  identity  satis- 
fied by  any  set  of  the  X-  and  Y-functions  [cf.  Chandrasekhar,  1950, 
p.  187;  Sekera,  1966b,  Equation  (54)].  If  we  use  for  the  modified 
moments  the  following  notation; 


y x^(t)  dt,  = y Y^(t)  t"  dt 


(III-41) 


(III-42) 


The  modified  moments  x and  y can  be  expressed  in  terms  of  the 
ordinary  moments,  defined  in  Equation  (11-46)  by  substituting  in 
Equation  (III-41)  the  expressions  for  the  corresponding  functions.  In 
this  way  we  obtain 


for  : 


’'n  ■ 8 * ^“l^n+2  " ^^“l^n+4^’  ^ J = 


(III-43) 


J 


for  X^: 


(2) 


dt  = 0.70 


1 

" T6  ^^'^2^n+2  ^^2^n+4^’  ^ f " 

0 

(III-44) 


1 

for  X^:  ^n  " I ‘ ^‘^3^n+2^’  ^ ^ 

0 


(111-45) 

1 

for  X^:  x^  = I [(a^)^  - 2 J dt  = 0.50 

0 


(111-46) 


The  expressioi’.s  for  y^  are  obtained  from  Equations  (111-43)  through 
(111-46)  by  replacing  a by  B . Table  D gives  the  values  of  the 
left  side  of  Equation  (III-42)  correspondi ng  to  various  optical  thick- 
nesses for  the  functions  (i  =1,2,  4).  For  i = 3,  the  linear 

constraints  on  X3  and  lead  to  the  relationship 

Xq  ^ yo  = 1 (111-47) 


70 


The  values  of  Xq  + for  various  optical  thicknesses  are  also  given 
in  Table  D.  The  relations  in  Equations  (III-38),  (III-39),  (III-42), 
and  (III-47)  allow  accuracy  checks  of  the  moments  of  even  order  up  to 


p = 4 for  , X2,  '^2’  order  0,  1,2,  for  X^,  Y^, 

X^,  Y^  . To  check  the  accuracy  of  higher  moments,  one  can  use  the  fol- 
lowing relations  [cf.  Chandrasekhar,  1950,  p.  188;  Sekera,  1966b,  Eqs. 
(88)  and  (89)]; 


(1  - 


-I 


C;(t)t‘^dt  (III-48) 


for  , Y^  (i  = 1 , 2 , 4) , and 


(x£  y2)(l  " Xq  '*■  Yq)  (x-j  - y-|)  " ^ 


(III-49) 


for  X^,  Y^  . 
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COMPARISON  OF  THE  ACCURACY  OF  THE  K.-  AND  L. -FUNCTIONS  FOR  u = 1.00 


IN  VARIOUS  TABLES  WITH  THE  USE  OF  THE  RELATION  IN  EQ.  (III-30) 


T 

Ki  - K3 

BHa 

>-2  ■ ^4 

Authors 

0.15 

-.00011 

+.00005 

-.00010 

-.00017 

Chandrasekhar  and  Elbert, 
1954 

.00001 

.00000 

.00001 

-.00000 

Sekera  and  Blanch,  1952 

.00000 

.00000 

.00000 

.00000 

Present  tables 

0.25 

-.00011 

+.00009 

-.00009 

-.00016 

Chandrasekhar  and  Elbert, 
1954 

-.00006 

-.00000 

-.00005 

+.00000 

Sekera  and  Blanch,  1952 

.00000 

.00000 

.00000 

.00000 

Present  tables 

1.00 

-.00967 

-.00018 

-.00535 

+.00179 

Chandrasekhar  and  Elbert, 
1954 

-.00054 

+.00001 

-.00025 

-.00016 

Sekera  and  Blanch,  1952 

.00000 

.00000 

.00000 

.00000 

Present  tables 

I 
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Table  B 

ACCURACY  CHECK  OF  THE  MOMENTS  OF  THE  K.-  AND  L. -FUNCTIONS  BY  USE  OF 
THE  RELATIONS  IN  EQS.  (111-32),  (III-33),  (III-34),  AND  (III-35) 


T 

M^ 

"3 

^^4 

0.15 

1.00000 

1.00000 

0.00000 

-0.00000 

0.50 

1.00001 

1.00000 

0.00000 

0.00001 

1.00 

1 . 00000 

0.99999 

-0.00000 

0.00000 

4.00 

1.00000 

1.00000 

-0.00002 

-0.00005 

16.00 

1.00001 

0.99998 

0.00005 

0.00030 

100.00 

0.99999 

0.99999 

-0.0010 

-0.0015 
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Table  C 

ACCURACY  CHECK  OF  THE  MOMENTS  OF  THE  FUNCTIONS 
X.  AND  Y.  (i  = 3,  4)  SATISFYING  THE  RELATIONS 
IN  EQS.  (III-38)  AND  (III-39) 


T 

For  i = 3: 
^1 

For  i ^ 4: 
^2 

0.15 

0.99999 

0.99999 

0.25 

1.00000 

1.00001 

1.00 

0.99999 

0.99997 

4.00 

1.00007 

0.99992 

16.00 

1.00004 

1.00024 

100.00 

1.0000 

1.0000 

Table  D 

ACCURACY  CHECK  OF  THE  MODIFIED  MOMENTS  OF  THE  FUNCTIONS  X.,  Y. 

(i  = 1,  2,  4)  FROM  THE  RELATION  IN  EQ.  (III-42)  AND  OF  THE  FUNCTIONS 
X3  AND  Y3  FROM  THE  RELATION  IN  EO.  (III-47) 


2^0  - ^0  ^ ^0  ''O  ^ ^0 


T 

i = 1 

i = 2 

i = 4 

i = 3 

0.15 

0.69998 

0.70001 

0.49999 

0.99998 

0.50 

0.69998 

0.70001 

0.50001 

0.99997 

1.00 

0.70001 

0.70001 

0.49998 

1 .00003 

2.00 

0.70000 

0.70000 

0.50001 

0.99999 

4.00 

0.70002 

0.70001 

0.49998 

0.99997 

8.00 

0.70001 

0.69999 

0.50002 

1.00000 

16.00 

0.70001 

0.69999 

0.50002 

0.99999 

100.00 

0.70001 

0.70000 

0.50002 

0.99996 

0.700000 


0.700000 


0.500000 


1 .000000 


IV. 


SOME  APPLICATIOnS  OF  THE  TABLES 


f 

The  purpose  of  this  chapter  is  to  briefly  illustrate  the  type 
of  infonnation  v;hich  can  be  obtained  from  the  solution  to  the  radiative 
transfer  problem  as  represented  by  the  tables  of  X-  and  Y-  and  K-  and 
L-functions.  A more  complete  description  can  be  found  in  the  published 
papers  of  Kahle  (1968a,  1968b)- 

A . INTENSITY  OF  EMERGENT  RADIATION 

In  this  section  we  shall  examine  the  characteristi cs  of  the  in- 
tensity of  the  radiation  emerging  from  the  upper  and  lower  boundaries 
of  a plane-parallel  atmosphere  for  a range  of  external  parameters  (solar 
elevations,  ground  reflectivities,  angle  of  emergence)  for  the  entire 
domain  of  the  optical  thickness.  Chandrasekhar  and  Elbert  (1954)  and 
also  Coulson  et  al.  (1960),  published  tables  of  the  intensity  and  po- 
larization parameters  (Stokes  parameters)  for  radiation  emerging  from 
the  bottom  and  top  of  the  atmosphere  for  i 1 1-  Sekera  (1957)  has 
illustrated  graphically  many  of  the  features  of  the  polarization  of  the 
downward  radiation,  and  Coulson  (1959)  has  shown  graphs  of  the  upward 
intensity  and  polarization  for  the  same  optical  depths.  Dave  and 
Furukawj  (1966)  have  also  shown  a few  of  the  main  features  of  the  in- 
'iri.it>  jf'J  polarization  values  graphically  for  larger  optical  thick- 
1).  ver,  the  results  presented  here  are  the  first 

't  th('  intensity  c.f  radiation  emerging 


from  a plane-paral  )e1  Rayleigli-scattoring  atinosphere  of  large  optical 
thi ckncss . The  optical  thickness  of  a Rayleigh  atmosphere  is  strongly 
dependent  upon  the  wavelength,  ■ [Deinnendjian,  1955].  The  intensity 
depends  upon  the  optical  thickness,  t (and  thus  indirectly,  the  wave- 
length, A ),  upon  the  angles  of  incident  radiation,  (zenith  angle) 
and  (1q  (azimuth  angle),  upon  the  angles  of  emerging  radiation,  9,  p 
and  upon  the  ground  reflectivity,  A.  We  assume  Lambert  reflection 
and  allow  the  reflectivity  of  the  ground  to  range  from  A = 0 (total 
absorption)  to  A = 1 (total  reflection).  The  effect  of  ground  re- 
flectivity is  to  add  a term  to  the  intensity,  I , which  consequently 
changes  the  percent  of  polarization,  fO  , but  the  polarization 
parameters  Q and  U are  unchanged,  since  a depolarizing  reflection 
is  assumed  (cf.  Chapter  III,  Section  C).  Only  the  intensity,  I , is 
examined  in  this  section. 

An  incident  radiation  of  unit  intensity  per  unit  area  is 
assumed  throughout  this  study.  Figures  1 and  2 illustrate  the  plane- 
parallel  approximation.  The  intensity  vectors  of  the  scattered  radia- 
tion emerging  from  the  upper  and  lower  boundaries  of  the  atmosphere  for 
a specific  set  of  parameters,  t = 1.0,  Oq  = 53.13°  (pq  h cos  0q  = 
0.6),  A = 0,  are  shown  in  the  principal  plane--the  plane  containing  the 
sun,  the  point  of  observation,  and  the  local  vertical.  The  direct 
solar  radiation,  which  decreases  exponentially  with  optical  thickness 
as  the  radiation  travels  through  the  atmosphere,  is  not  included.  Fig- 
ures la  and  2a  show  diagramatical ly  how  the  directions  of  the  upward 


! 
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and  downward  emergent  radiation  are  related  to  the  model.  Figures  lb 
and  2b  show  the  same  data  in  a form  suitable  for  comparison  and  inter- 
pretation. The  variation  with  ^ , the  azimuthal  angle,  will  also  be 
considered . 

The  height  of  the  atmospheric  slab  in  these  figures  is  merely 
diagramatic.  No  geometric  height  is  determined  since  the  radiation  is 
a function  of  optical  thickness  only.  For  downward  radiation  the  actual 
height  is  practically  immaterial  except  for  the  effect  it  has  on  the 
validity  of  the  plane-parallel  approximation.  When  making  measurements 
of  upward  radiation,  however,  one  must  be  high  enough  above  the  planet 
to  be  above  substantially  all  of  the  atmosphere,  yet  at  the  same  time 
low  enough  that  the  planetary  surface  below  is  well  approximated  by  a 
plane  surface. 

Alternatively,  for  observations  from  great  distances,  the  emerg- 
ing radiation  can  be  approximated  by  parallel  radiation.  Sekera  and 
Viezee  (1964)  made  calculations  of  the  intensity  and  polarization  of 
such  radiation  for  1,  still  using  the  plane-paral  lei  approximation 
but  moving  the  point  of  tangency  to  each  point  of  the  observed  disk. 

This  method  becomes  less  reliable  as  the  ratio  of  the  geometric  thick- 
ness of  the  atmosphere  to  the  planetary  radius  increases.  For  inter- 
mediate distances,  where  neither  approximation  is  valid,  as  for  instance 
at  the  altitude  of  the  ATS  satellites  at  three  earth  radii,  each  case 
must  be  calculated  with  its  own  particular  geometry. 

We  can  investigate  how  the  intensity,  I , varies  with  optical 
thickness,  direction  of  incident  light,  and  ground  reflectivity.  We 
will  consider  hero  the  downward  case,  the  intensity  of  the  transmitted 
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skylight  as  would  be  seen  by  an  observer  at  the  planet's  surface, 
assuming  unit  incident  radiation  per  unit  area,  at  the  top  of  the 
atmosphere. 

With  the  sun  in  the  zenith  the  intensity  must,  by  definition, 
be  symmetric  around  the  local  vertical;  there  is  no  azimuthal  (;) 
variation.  Therefore,  only  half  of  the  principal  plane  need  be  illus- 
trated. Figure  3 shows  the  variation  with  both  optical  thickness,  t , 
and  ground  reflecti vi ty , A,  for  this  case  of  the  sun  in  the  zenith. 
The  left  half  of  the  diagram  presents  the  intensities  for  A = 0 and  the 
1 values  shown;  the  right  half  shows  the  intensities  for  the  san’e  : 
values  and  A = 0.8.  In  both  cases  we  see  the  maximum  intensity  in  the 
zenith  is  at  about  t = 2,  while  at  the  horizon  the  maximum  is  at  about 
T = 1 . For  small  i the  intensity  is  greatest  at  the  horizon,  but 
by  T = 1 the  intensity  becomes  greatest  at  the  zenith.  Even  for 
T = 100  there  is  still  a significant  amount  of  radiation  that  pene- 
trates through  the  atmosphere  to  the  planet's  surface.  The  effect  of 
a high  ground  reflectivity  is  to  increase  the  intensity  substantially 
(by  roughly  a factor  of  2)  for  all  values  of  i , and  at  all  zenith 
angles.  For  moderate  optical  thickness  (x  = 0.5  to  x = 4)  this 
increase  of  intensity  is  considerably  greater  near  the  horizon  than 
near  the  zenith,  significantly  changing  the  shape  of  the  intensity 
curves . 

When  the  sun  is  not  in  the  zenith,  the  axial  symmetry  about  the 
zenith  is,  of  course,  lost,  the  only  obvious  symmetry  remaining  being 
across  tiie  principal  plane.  As  an  example  of  moderate  solar  zenith 
angle,  v/e  will  examine  the  radiation  when  0^  = 53.13"  (i.q  = 0.6). 
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The  inLensiLy  in  Ihe  principal  plane  with  i as  a parameter  is  shown 
in  Figure  4a  for  A = 0,  and  in  Figure  4b  for  A - 0.8.  Both  sets  of 
curves  show  the  same  features,  differing  more  in  magnitude  than  in 
shape,  with  the  high  ground  reflectivity  again  almost  doubling  the  in- 
tensity. For  low  optical  thicl.ness,  tite  intensity  is  greatest  near  the 
horizon.  For  increasing  optical  tiiickness  these  two  peaks  of  intensity 
move  toward  the  zenith  until,  somewhere  between  : = 2.0  and  i = 4.0, 
the  two  maxima  become  one,  slightly  to  the  sunward  side  of  the  zenith. 

By  1 = 8,  all  asymmetry  from  one  side  of  the  zenith  to  the  otiier  is 
essential  1>  gone,  with  the  peak  intensity  in  the  zetitli.  Coriparing 
Figures  3 and  4,  we  see  that  the  horizon  brightening  at  low  ’ is  more 
pronounced  wiien  the  sun  is  not  in  the  zenith,  while  the  average  total 
intensity  is  naturally  lower. 

Figures  5a  and  5b  stiuw  the  azimuthal  variation  ’ntens'*^.  ♦ n 
two  of  tiie  cases  already  cxai-i ne J . Ttiese  figures  a proje.t  ;:>■  f 

tiie  sky  onto  a pKme:  tiie  diar  eter  of  ttie  sei;  i r i rcl;-.  are  of  en -rl 
zenit’n  distariLO,  and  tiie  radijl  lines  are  of  equal  u.'’  r.itb.  Orly  ’ '■* 
ttic  circle  is  shown  sitice  the  ra'liation  is  sy  Htrwal  at:u'  the  p'lr  , - 
pal  plane.  Tiie  contoni  . arc  lines  of  epi-il  I'  tensi'y.  Figu-a  , f;,, 

1 - 1.0,  Py  0.6  and  A ' 0 (cf  figure  4<i ) . show-  ho.-,  the  iiitin.it.. 
varies  witli  azimuth  \/hon  the  pcai  intetisity  is  quite  far  from  th* 
horizon,  figure  5h , for  * 0.15,  - 0.6  and  A 0.8  (cf  figure  4;) 

shows  the  otiiet  ex.tri'i.o.  wtinin  tlie  intensity  is  groMest  on  the  hori- 
izun.  [iutli  figuri's  siiov;  a lower  v.i  r i a t i(iri  of  ifiten.sity  with 
aziimit.h  angle  tiwni  with  zenith  angle,  liie  peak  of  intensity  in  figure 
5a  is  somewtiat  broader  than  it  is  high.  Tiiis  effect  is  much  more  pro- 
nounced in  Figure  5b.  The  veiy  liright  horizon  continues  all  tiic  way 
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around  witli  only  a relatively  slight  decrease  vnth  increasing  distance 
from  the  principal  plane. 

There  is  an  unexpected  feature  in  Figures  5a  and  5b  that  con- 
tinues through  all  the  calculations  of  botii  intensity  and  polarization. 
This  is  a symmetry,  on  the  horizon  only,  around  the  ; = 90°,  ; = 270° 
line,  that  is,  the  line  perpendicular  to  the  principal  plane,  through 
the  point  of  observation.  In  the  other  figures  one  can  also  see  this 
effect,  in  that  both  ends  points  of  each  curve  (■:•  = 0°  and  ; = 180°  , 
0 = 90°)  are  at  the  same  value.  The  reason  for  this  sym.T,etry  is  not 
clear.  Since  it  only  occurs  on  the  horizon,  where  the  plane-parallel 
approxiPiUti  on  loses  its  validity  for  a planetary  atir.osphere , it  will 
probably  not  be  observed  in  a real  atmosphere. 

The  results  for  only  a few  representative  values  of  the  solar 
zenith  angle  and  ground  reflectivity  have  been  illustrated  here,  bat 
the  intensity  of  the  eimargerit  radiation  can  bo  dete*'mi :;c J from  the 
tables  in  the  Appendix  for  any  desired  values  of  these  parameters.  Simi- 
lar results  for  the  upward  intensity,  as  would  be  measured  by  a lew- 
orbiting  Satellite,  can  also  be  derived  from  the  tct.les,  but  will  not 
be  illustrated  here.  examples  of  the  ui-war-c;  intensity  are  given 

in  the  paper  by  Kaiile  (19CCb),  which  is  reproduced,  here-  in  Appendix  B. 

A Raylcigh-scattcring  at;  osphere  is  only  an  approx ir.ati on , of 
course,  to  any  real  planetary  atr.ospiiore . Ttie  a.ppl icabi  1 i ly  of  tf.is 

approximation  and  its  1 in'i  t atic>ns  have  been  tiiscussed  in  detail  else- 
where (sec,  for  example,  Ctiandrasekhar , 1950;  Deirmondji an , 1 957,  1959; 
Rozenborg,  19G6;  Sebera,  1957  ).  It  is  h(;|ud  that  even  when  a rc.al  at- 
mosphere differs  considerably  from  a Rayloigii-scat  tering  atmosi'liore , 
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useful  inforiiidtion  can  be  gainod  by  comparing  it  with,  and  noting  the 
difference  from,  the  theoretical  model. 


B . GLOBAL  RAD I A1 ION 

Global  radiation,  in  the  sense  used  here,  is  the  total  radia- 
tion of  a given  wavelength  reaching  a unit  surface  area  of  a planet, 
and  includes  both  the  directly  transmitted  solar  radiation  and  the  dif- 
fusely scattered  radiation  from  the  entire  sky.  Tiie  diffuse  part  is 
just  the  hemispherical  integral  of  the  directional  intensity  illus- 
trated in  the  previous  section.  The  relative  global  radiation  is 
usually  expressed  as  a fraction  o;  the  radiation  incident  at  the  top 
of  the  atir.ospliere  within  a given  spectral  range.  Satellites  now  can 
Pleasure  the  upward  global  radiation  emerging  from  a unit  surface  area 
at  the  top  of  the  atmosphere,  so  v.'C  distinguish  between  upward  and 
dowruvard  global  radiation. 

We  assui.e  an  incident  flux  of  parallel  radiation  frci  tiie  sun  of 
T!  units  [)cr  unit  area  perpendicular  to  the  direction  of  propagation; 
i.e.,  Fq  = 1 . A unit  area  of  the  top  of  the  atmosphere  thus  receives 
Wg-'  units,  wliorc  ..g  is  the  cosine  of  the  solar  ;’oni  th-di  stance  -g. 
The  global  I'ddiation  can  be  found  from  Lquations  (III-?3)  through 
(111-26). 

In  iigurc  6 wo  show  on  a log-log  scale  the  sun  radiation  S,  the 
downwari.1  diffuse  sly  radiation  11^  , and  the  upward  diffuse  I'adiation 
11^  as  a function  of  optical  thickness  : for  vettical  incidence  of 

solar  radiation  and  no  ground  ref l('Cti  vi  ty . The  sum  of  these  three 
components  for  any  givc'ti  optical  thickness  must  by  definition  equal  the 
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incident  radiation  > of  in  this  case  just  ^ . For  the  continua- 

tion of  these  curves  in  the  direction  of  smaller  ; , see  Deirmenjian 
and  Sekera  (1954).  It  can  be  seen  that  when  the  optical  thickness 
reaches  4,  the  direct  solar  radiation  is  an  order  of  magnitude  smaller 
than  the  downward  diffuse  radiation,  and  with  increasing  optical  thick- 
ness it  rapidly  becomes  negligible.  Although  the  upward  radiation  is 
approaching  ^ asymptotically,  there  is  still  a measurable  amount  of 
downward  diffuse  radiation  (almost  2 percent)  at  ■ = 100. 

As  with  the  intensity  in  the  previous  section  we  will  examine 
only  a couple  of  examples  of  how  the  global  radiation  varies  with  the 
parameters,  solar  zenith  angle  and  ground  reflectivity.  More  examples 
of  these  variations  are  given  in  the  paper  by  Kahle  (1968a),  and  again, 
of  course,  calculation  over  the  complete  range  of  parameters  is  possible 
from  the  functions  given  in  the  tables  of  the  Appendix. 

The  variation  of  the  diffuse  downward  radiation  H^,  for  an 
angle  of  incidence  Oq  and  no  ground  reflection,  is  shown  in  Figure  7. 
This  radiation  shows  an  interesting  change  with  solar  zenith-distance. 
The  peak  flux,  at  about  i - 2 for  vertical  incidence,  moves  to 
smaller  a with  increasing  solar  zenith-distance,  until  for  Oq  ■-  88° 
it  has  moved  below  i = 0.15  . The  peak  is  caused  by  the  balance  be- 
tween the  increased  scattering  as  path- length  through  the  atmosphere 
increases,  and  the  reduction  of  the  total  light  penetrating  to  lower 
levels  as  optical  thickness  increases. 

We  next  consider  the  effect  of  ground  reflectivity,  again 
assuming  Lambert's  law.  Figure  8 shows  the  downward  diffuse  radiation 
as  a function  of  t for  various  values  of  A and  for  Uq  = 1.0. 
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(Other  values  of  Uq  merely  reduce  the  flux  in  a way  similar  to  that 
shown  in  Figure  7. 

This  downward  radiation  shows  a somewhat  unexpected  feature; 
for  very  large  surface  reflectivities  there  is  more  radiation  impinging 
upon  the  planet's  surface  at  the  bottom  of  a thick  atmosphere  than  is 
incident  at  the  top  of  the  atmosphere.  For  A = 0.95  this  is  true  be- 
tween T = 2 and  = 7.  For  total  ground  reflectivity,  A = 1 , this 
phenomenon  occurs  for  all  larger  tfian  1.5  , and  for  i larger  than 
8 the  amount  of  downward  diffuse  radiation  becomes  completely  indepen- 
dent of  optical  thickness.  It  seems  paradoxical  at  first  that  the 
amount  of  radiation  incident  at  the  surface  of  a planet  that  is  covered 
by  an  atmosphere  so  thick  that  a negligible  amount  of  direct  solar  ra- 
diation can  penetrate  it  can  exceed  the  amount  of  radiation  incident  at 
the  top  of  the  atmosphere.  Hov;ever,  this  is  just  because  the  radiation 
which  does  reach  the  bottom  surface  becomes  partly  trapped  between  the 
perfect  reflector  below  and  the  diffusively  reflecting  atmosphere  above, 
undergoing  multiple  reflections  before  eventually  escaping  to  space. 

That  there  is  no  violation  of  energy  conservation  can  be 
demonstrated  by  considering  the  fluxes  involved  in  the  energy  balance. 
The  total  radiation  incident  at  the  lower  boundary  does  not  represent 
a flux  through  it  unless  A = 0. 

These  results  have  been  presented  in  a general  form  so  that  they 
can  be  applied  to  a variety  of  problems.  The  principal  parameter, 
optical  thickness,  is  related  to  the  wavelength  of  the  radiation  by 


6,  dz  (IV-1) 

z=0 

where  z is  the  vertical  coordinate.  The  extinction  coefficient, 

, is  proportional  to  Thus  one  could  equally  well  assume 

that  v/avelength  is  the  independent  variable  in  the  flux  curves.  For 

example,  the  optical  thickness  of  the  earth's  molecular  atmosphere 

0 

ranges  from  about  t = 0.454  at  A = 3750A  (blue)  through  • t = 0.15 
0 0 
for  A = 4950A  (green)  to  x = 0.0173  at  A = 8350A  (red) 

[Dei rmendjian , 1955].  In  Figure  8 v/e  see  that  the  downward  diffuse 
flux  at  X = 0.5  (blue)  is  several  times  larger  than  at  x = 0.15 
(green),  and  from  the  work  of  Dei rmendjian  and  Sekera  (1954)  we  see 
that  the  flux  continues  to  decrease  down  into  the  red.  Thus  most  of 
the  energy  of  the  skylight  will  be  in  the  blue  region  of  the  spectrum. 
(Deirmendjian  and  Sekera  obtained  quantitative  results  by  including  the 
spectral  distribution  of  the  incident  solar  radiation.)  If  one  were 
next  to  assume  a much  thicker  atmosphere,  several  times  as  thick,  op- 
tically, as  the  earth’s,  then  each  region  of  the  spectrum  would  shift 
to  a correspondingly  higher  optical  thickness.  This  would  place  the 
visible  part  of  the  radiation  beyond  the  peak  of  the  flux  curve  in 
Figure  7.  Hence  the  spectral  characteristics  of  the  diffuse  sky  radia- 
tion would  bo  changed,  and  a much  greater  portion  of  the  emerging 
energy  would  lie  in  the  red. 
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To  make  accurate  predictions  regarding  a real  planetary  atmos- 
phere on  the  basis  of  the  foregoing  analysis,  it  is  necessary  to  know 
how  the  atmosphere  departs  from  idealized  conditions.  One  must  con- 
sider how  closely  the  atmosphere  corresponds  to  a Raylei gh-scattering 
atmosphere,  and  how  the  reflection  characteri sti cs  of  its  lower  surface 
(ground  or  cloud  layer)  depart  from  Lambert's  law.  Conversely,  as  more 
data  are  gathered  on  the  optical  properties  of  an  atmosphere,  the  dis- 
crepancy between  the  observed  and  the  predicted  values  will  provide  in- 
formation about  other  physical  properties,  such  as  atmospheric  turbid- 
ity, true  absorption  and  re-emission  of  light,  variations  of  the  com- 
position of  the  atmosphere  with  height,  and  the  nature  of  the  ground 
surface. 
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APPENDIX 

TABLES  or  SCATTERING  FUNCTIONS 


X-  and  Y-Functions 

Moments  of  the  X-  and  Y-Functions 

K-  and  L-Functions 

Moments  of  the  K-  and  L-Functions 
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Table  1 v-  MiD  Y-FUMCTIONS  (t  = 0.15) 


IJ- 

^1 

^2 

^2 

0.00 

0.  lOOOOF 

01 

0. 

0. lOOOOE 

01 

0. 

0.0^ 

0.  1032Ht 

01 

0. t6773t- 

-0  1 

0. 1C209E 

01 
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-01 

0.0'. 

0. 10662( 

01 

0. 66698F  ■ 

-01 

0.  10361E 

01 
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-01 
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0. 1071  7fc 

01 

0.  1 302U 

00 

0.  10677E 

01 

0. 1 1 731E 

00 

0.08 

0.  10H62F 

01 

0.216821 

00 

0.  10666E 

01 

0. 19777E 

00 

0.  10 

0.  10938F. 

01 

0.296611 

00 

0.  10636E 

01 

0.27513E 

00 

0.  12 

0. 1 1012F 

01 

0.367961 

00 

0.  1C687F 

01 

0. 36663E 

00 

0.1'. 

0. 1 10  ?2F 

01 
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00 
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01 
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00 

0.  16 

0 . 1 1 1 2 If 

01 
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00 

0.  10766E 

01 

0.66902E 

00 

0.  16 

0.  1 1 16U 

01 

0, 63696F 

00 

0.  10  796E 

01 

0. 60667E 

00 

0.20 
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01 
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00 
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0.66638E 

00 

0.22 
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0. 780188 

00 

0. 1 22968 

01 

0.  3 3 1 698- 

00 

0.66 

0.  186671 

01 

0.806 1 1 8 

OO 

0.  123168 

01 

0. 366288 

00 

0.68 

0. 186718 

01 

0.827618 

00 

0.  12  3188 

01 

0. 168668 

00 

0.  10 

0.  187  708 

01 

0.860  168 

00 

0.12  1688 

01 

0.171 768 

00 

0.  7? 

0. 188668 

01 

0.872668 

00 

0.  12  3708 

01 

0. 366638 

00 

0.  76 

0. 189698 

01 

0. 896618 

00 

0. 123978 

01 

0.  197268 

00 

0.  76 

0.  190681 

01 

0.916698 

00 

0. 126168 

01 

0.609618 

00 

0.  78 

0.191361 

01 

0.936668 

00 

0.  1263  18 

01 

0.621 768 

00 

0.80 

0. 19?1  7f 

01 

0.966  718 

00 

0. 126618 

01 

0.613668 

00 

0.8? 

0.  19?9'I8 

01 

0.976698 

00 

0.  126678 

01 

0.666298 

00 

0.86 

0.193  r.t 

01 

0.996  798 

00 

0.  1 268  18 

01 

0. 666698 

00 

0.86 

0. 19661 L 

01 

0.101668 

0 1 

0.  1 26988 

0 1 

0.68.  7868 

CO 

0.8  8 

0.  196?lt 

01 

0.  101 108 

01 

0.  I 26 1 18 

01 

0. 6 78  798 

00 

0.  <J0 

0. 196961 

01 

0. 106108 

0 1 

0.  126278 

01 

0.689608 

00 

0.9? 

0. 196671 

01 

0. 106868 

01 

0.126618 

01 

0.6999 98 

on 

0.96 

0.  197?1U 

01 

0.  1086  78 

0 1 

0 . 1 2 6 6 S 8 

01 

0.610268 

00 

0.96 

0.  19  77?l 

01 

0.  1 1 0268 

0 1 

0.  126688 

01 

0.620308 

00 

0.98 

0 . 1 9 8 6 •.  1 

01 

0. 1 1 1888 

01 

0.  126808 

01 

0. 630  1 38 

00 

1 .00 

0.199161 

01 

0.113678 

01 

0. 1 26928 

01 

0.639778 

00 

I 
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Table  1 (coni.)  (f  - 2,0) 


^2 

^2 

0.00 

0.  lOOOOf: 

01 

0. 

0. 100001 

01 

0. 

0.0? 

0. I0366t 

01 

0.128371 -C2 

0. 102311 

01 

0.  171911-02 

O.O'. 

0.1061 3f 

01 

0.268  321 -02 

0.106121 

01 

0.  366261  - 02 

0.  06 

0.  1 08  3St 

01 

0.619821-02 

0.106721 

01 

0.660191-02 

0.08 

0.  1 10331 

01 

0.682661-02 

0.  10  7 1 

01 

0.  766601  - 02 

0.  10 

0.  1 1216E 

01 

0.  766  (71-02 

0 . 10  8 6 71 

01 

0.976261-0? 

0.  i? 

0, 1 1 301 t 

01 

0.963631-02 

0. 1 09861 

01 

0.  1206  11-01 

0.  K 

0.  1 1 63  7f 

01 

0. 1 162  (1  -01 

0, 1 1 1091 

01 

0.  166631-01 

0.16 

0.  1 1 68  )E 

01 

0. 1 36661 -01 

0. 1 1 2261 

01 

0.  170191-01 

0.  IH 

0.118211 

01 

0. 168681  -0  1 

0.  1 1 5 5 71 

01 

0.197161-01 

0.  ?0 

0. 1 1 9621 

01 

0.183121-01 

0.  1 1 6631 

01 

0.226761-01 

0.22 

0.120761 

01 

0.209  761  -01 

0.  1 1 6661 

01 

0.266161-01 

n.  2'. 

0.  121931- 

01 

0.23H731 -01 

0.116661 

01 

0.28-3691-01 

0.26 

0.  123061 

01 

0.270391-01 

0.  1 1 7 381 

01 

0.323611-01 

0.28 

0.1261 31 

01 

0.  306061--0  1 

0.118301 

01 

0.361211-01 

0.  30 

0.  126161 

01 

0. 3630V1 -01 

0. 1 19141 

01 

0.601811-01 

0.  32 

0.126161 

01 

0.  3866  7f -0 1 

0. 120021 

01 

0.666621 -01 

0. 3<f 

0.  12  708f 

01 

0.629931 -01 

0. 120861 

01 

0.697691-01 

0.  36 

0. 127991 

01 

0.679071-01 

0. 121661 

U 1 

0. 66  3231-01 

0.  38 

0. 128861 

01 

0.6322  31-0  1 

0. 122601 

U1 

0.697341-01 

0.60 

0.129  701 

01 

0.6892  71  -01 

0.123161 

01 

0. 666001 -01 

0.62 

0. 1 30611 

01 

0.66021 101 

0. 123861 

01 

0.716121-01 

0.66 

0. 131291 

01 

0. 7 16961-0  1 

0. 126661 

01 

0.780631-01 

0.66 

0.  1 32061 

01 

0.7832  71-0 1 

0.  126261 

01 

0.8683  71-01 

0.  68 

0.  1 32  761 

01 

0.  866 1 6F-0 1 

((.  126891 

01 

0.919291-01 

0.60 

0- 1 33661 

01 

0.930201-01 

0. 126621 

01 

0. 993061-01 

0.  62 

0. 136161 

01 

0.100861  00 

0.127161 

01 

0.106971  00 

0.66 

0.  1 36  791 

0 1 

0.108961  00 

0.  1 2 7761 

01 

0.116881  no 

0.  66 

0.  i 36621 

0! 

0.  117  111  00 

0.  1 28  371 

01 

0.123061  00 

C.  68 

0.  136031 

01 

0.126961  00 

0.  128881 

0 1 

0.131601  00 

0.60 

0.  1 36631 

01 

0.136781  00 

0. 129631 

01 

0.139971  on 

0.62 

0. 137201 

01 

0.163811  00 

0. 129961 

01 

0.168711  00 

0.66 

0.  1 37  761 

01 

0.163021  00 

0. 130671 

01 

0.167611  00 

0.66 

0.  1 38291 

01 

0.162601  00 

0. 1 30981 

01 

0.166661  00 

0.68 

0.138821 

01 

0.171901  00 

0.131661 

01 

0.  1 763  1 1 (10 

0.  70 

0.  I 39  321 

01 

0.181661  00 

0.131 961 

01 

0.1B6071  00 

0.  72 

0. 1 39821 

01 

0.191271  00 

0. 132601 

01 

0.196631  00 

0. 76 

0. 160291 

01 

0.201091  00 

0.132861 

01 

0.20  3 8 61  00 

0.  76 

0.  160  761 

01 

0.210(81  00 

0.  1 3 3291 

01 

0.213361  00 

0.  78 

0.161211 

01 

0.220931  00 

0.133711 

0 1 

0.227881  00 

0.80 

0.161 661 

01 

0.230931  00 

0.136131 

0 1 

0.232661  00 

0.82 

0.  1620  71 

01 

0.260961  00 

0. 1 36631 

01 

0.262061  00 

0.  86 

0.  162691 

01 

0.261001  00 

0.  1 369  31 

01 

0.261661  00 

0.86 

0. 162891 

01 

0.261(161  Oil 

0.136311 

01 

0.261261  on 

0. 88 

0.  163281 

01 

0.2  7 1 101  00 

0. 1 36691 

0 1 

0.270861  00 

n.<70 

0.163671 

01 

0.281161  00 

0. 1 36061 

01 

0.780661  00 

0.92 

0.  166061 

01 

0.2911 71  no 

0. 1 36611 

01 

0.290021  CO 

0.96 

0.  1 66601 

01 

0.  301  1 71  00 

0. 1 36761 

01 

0.299661  00 

0.96 

0.  166  761 

01 

0.311161  00 

0.  1 3 7101 

0 1 

0.3090  71  00 

0.98 

0. 166101 

01 

0.  3210  71  00 

0.  1 3 76  31 

01 

0.  318631  00 

i .00 

0. 1 66661 

01 

0.330961  00 

0.  1 3 7 761 

01 

0.327961  00 
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T ah  1 1 

1 (coat.) 

(T 

2.0) 

S 

"3 

1 , 
•4 

0.00 

0. lOOOOl 

01 

0. 

0.  lOOOOl 

C 1 

0. 

0. 0^ 

0. 107661 

0 1 

0.696601- 

-02 

0. 102921 

ni 

0.402741-03 

0. 0<. 

0.11  1641 

01 

0. 147441 

-01 

0.104881 

0 1 

0.643091-03 

0.06 

0 . 1 1 8 b U 

01 

0.233061 

-0  1 

0.  1064  71 

01 

0.  1 32  2 31-02 

0.08 

0.  1?371F 

01 

0.  326  7 1 t- 

-01 

0.107881 

0 1 

0.  184061-02 

0.10 

0.128331 

01 

0 . 4 2 6 2 8 1 

-0  1 

0.  10  7 1 IF 

U 1 

0.2401 11-02 

0.  1? 

0.  132  731 

01 

0.6316  31 

-01 

0. 1 10221 

01 

0.  500  791-02 

0.  u 

0. 1 36961 

01 

0. 644  6 11 

-01 

0.  1 1 ! 2 31 

01 

0.  3664  71-02 

0.16 

0.141021 

01 

0. 764161- 

-01 

0.112161 

01 

0.436221-02 

0.  18 

0. 144961 

01 

0.890201- 

-0  1 

0.113011 

01 

0.616931-02 

0.20 

0.  1 48  761 

01 

0.  102  311 

00 

0.113811 

01 

0.6C606E-02 

0.22 

0.  1624  31 

01 

0. 1 16271 

00 

0.11 4661 

01 

0. 7046  71- 02 

0.2^ 

0. 166011 

)1 

0. 1 30931 

00 

0.  1 1 6261 

01 

0.019401-02 

0.2fj 

0. 1 69481 

01 

0. 146281 

00 

0.116911 

01 

0.962611-02 

0.26 

0. 162861 

01 

0.  162  341 

00 

0.  1 166  31 

0 1 

0.  1 10761-01 

0.  iO 

0. 166141 

01 

0. 1 791 11 

00 

n.  1 1 7121 

01 

0.128811-01 

0.  32 

0. 169341 

01 

0.  19667,1- 

00 

0.  11  7681 

01 

0.  149641-01 

0.  34 

0. 172461 

01 

0.214711 

00 

0. 118201 

01 

0.173641-01 

0.36 

0.176481 

01 

0.233491 

00 

0. 118701 

01 

0.200791-01 

0.  38 

0.  1 78431 

01 

0.262921 

00 

0.119181 

01 

0.231411-01 

0.  40 

0.101 301 

01 

0.272911 

00 

0.11  7641 

01 

0.266411-01 

0.42 

0. 184091 

01 

0.293461 

00 

0. 120071 

01 

0. 302881-01 

0.44 

0. 1 86811 

01 

0. 314611 

00 

0. 1 20491 

01 

0.  343761-01 

0.46 

0. 109461 

01 

0.336011 

00 

n. 1 20697 

01 

0.  38  7 961  -01 

0.48 

0.  t'72041 

01 

0. 357921 

00 

0.  12  12  71 

01 

0.4  364  3f-0I 

0.60 

0.194661 

01 

0.380181 

00 

0. 121641 

01 

0.466911-01 

0.62 

0. 197001 

01 

0.402771 

00 

0.121  971 

01 

0.639461-01 

0.64 

0.  199381 

01 

0.426621 

00 

0.122331 

01 

0. 696681- 01 

0.66 

0. 201 701 

01 

0.448701 

00 

0. 122661 

01 

0.664671-01 

0.  68 

0. 203961 

01 

0.471961 

00 

0. 122961 

01 

0. 716841-01 

0.60 

0.206161 

01 

0.496  361 

00 

0.  12  3261 

01 

0.  779361-01 

0.  62 

0. 208  31 1 

Cl 

0. 61H001 

00 

0.  1 2 3661 

01 

0.844961-01 

0.64 

0.210401 

01 

0.642461 

00 

0.  1 2 3841 

0 1 

0.912391-01 

0.  66 

0.212441 

01 

0.666091 

no 

0.124111 

01 

0.931671-01 

0.  68 

0.  21 4431 

01 

0. 689731 

CO 

0.124371 

01 

0.106221  GO 

0.  70 

0.216361 

01 

0.613361 

00 

0. 124621 

01 

0.112431  00 

0.72 

0.218261 

01 

0.636931 

00 

0.  1248  71 

01 

0.119761  00 

0.  74 

0.220101 

01 

0.660461 

00 

n. 126101 

01 

0.127201  OO 

0.  76 

0.221891 

01 

0.68  3881 

00 

0.  1 263  31 

01 

0.  1 34  731  00 

0.  78 

0.22  3661 

01 

0. 707201 

00 

0. 126661 

01 

0.142341  00 

0.80 

0. 226361 

01 

0.730411 

00 

0.126771 

01 

0.160011  00 

0.8? 

0.227031 

01 

0.  763471 

00 

0.126981 

01 

0.167761  00 

0.84 

0. 228661 

01 

0.  776  391 

00 

0 . 1 2 6 1 8 1 

01 

0.166631  00 

0.86 

0.2  30261 

01 

0.799141 

00 

0.126371 

01 

0.173341  00 

0.88 

0.  231801 

01 

0.821711 

00 

0. 126671 

01 

0.  18  1 I'll  00 

0.90 

0.  2 33  321 

01 

0.844111 

00 

0. 126761 

01 

0.189061  UO 

0.9? 

0.  2 34801 

01 

0.866311 

on 

0.126931 

01 

0.196911  00 

0.94 

0.236261 

01 

0.888301 

00 

0.127101 

01 

0.204791  00 

0.96 

0.  2 3 7661 

01 

0.910101 

00 

0.  1 2 7 2 71 

0 1 

0.212661  CO 

0.98 

0. 2 39061 

01 

0.931671 

00 

0.127441 

01 

0.220611  oO 

1 . 00 

0.240401 

01 

o 

• 

c 

(T 

00 

0.12  7601 

01 

0. 228361  CO 
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(T  = 4.0) 


Table  1 (cent  .) 


^2 

^2 

0.00 

0.  100008 

01 

0. 

0. 100006 

01 

0. 

0.0? 

0. 10366E 

01 

0.  186736-03 

0.  102  326 

01 

0. 371686-03 

o.o<. 

0.  1C616E 

01 

0.  386006-03 

0.  1061 36 

01 

0. 767086-03 

0.06 

0.  108  36t 

01 

0.602616-03 

0.  105756 

01 

0.  1 1 8 706-02 

0.08 

0.  n 036t 

01 

0.836566-03 

0.  1072  36 

01 

0.  163066-02 

0.  10 

0.  112166 

01 

0. 10825C-02 

0.  108626 

01 

0.209786-0? 

0.12 

0.  1 1 386E 

01 

0.  1 36676-02 

0.  109926 

01 

0.259006-02 

0.  1<* 

0. 1 1661t 

01 

0.  162756-02 

0.111166 

01 

0.  310666-02 

0.16 

0.1168  n 

01 

0.  19261 6-02 

0. 112366 

01 

0.  365036-02 

0.  18 

0. 1 1 826t 

01 

0.226266-02 

0.  1 1 3676 

01 
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0,20 

0.  1 1 96 7t 

01 

0.257926-02 

0.  1 16556 

01 

0.681806-02 

0.22 

0.  12082E 

01 

0.293586-02 

0.  1 15596 

01 

0.566676-02 

0.2<. 

0. 122006 

01 

0.  331  666-02 

0.  1 16506 

01 

0.610216-02 

0.26 

0. 123166 

01 

0. 371686-02 

0. 1 1 7556 

01 

0.679126-02 

0.28 

0. 126226 

01 

0.616636-02 

0.  1 1 8686 

01 
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0.30 

0. 126266 

01 

0.659996-02 

0.119376 

01 

0.827266-02 

0.  32 

0. 126266 

01 

0.506511-02 

0.  120266 

01 

0.906856-02 

0.  36 

0.  127216 

01 

0. 560696-02 

0. 12 1096 

01 

0.990666-02 

0.  36 

0.  1281 36 

01 

0.616086-02 

0. 121906 

01 

0.107886-01 

0.  38 

0.  12901E 

01 

0.675926-02 

0. 122696 

01 

0.  117176-01 

0.60 

0.  129876 

01 

0.760306-02 

0. 123666 

01 

0.  126986-01 

0.62 

0.  130696 

01 

0.809836-02 

0. 126216 

01 

0. 1 37366-01 

0.66 

0. 1316 36 

01 

0.885066-02 

0.  126966 

01 

0.168366-01 

0.66 

0. 132266 

01 

0.966626-02 

0. 125666 

01 

0.  159976-01 

0.68 

0.133016 

01 

0.  105686-01 

0. 126336 

01 

0.  1 72  336-01 

0.60 

0.  1 33736 

01 

0.  1 15036-01 

0. 127006 

01 

0.  185626-01 

0.  62 

0. 1 36636 

01 

0.  125606-01 

0.  127656 

01 

0.  199366-01 

0.66 

0.136116 

01 

0.136616-01 

0.  128296 

01 

0.216106-01 

0.66 

0.  136766 

01 

0.168766-01 

0.128916 

01 

0.229766-01 

0.68 

0. 1 36606 

01 

0.  161816-01 

0. 129516 

01 

0.266356-01 

0.60 

0.137026 

01 

0.  175886-01 
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01 

0.263926-01 

0.62 

0.  137636 

01 

0.  191006-01 

0. 1 30686 

01 

0.282696-01 

0.66 

0.138216 

01 
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01 
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0.66 
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01 

0.226676-01 

0.131 796 

01 
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0.68 

0.  139366 

01 
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01 
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01 

0.262666-01 

0.132856 

01 
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01 

0.283166-01 
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01 
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0.76 
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01 
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01 
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01 
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01 
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01 
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01 
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0.655626-01 
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01 
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01 
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01 

0.  587651  -1)1 
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01 
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01 
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01 
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01 
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01 
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01 
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